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Abstract 


The Krichevskii parameter determines the sign and magnitude of the near-critical divergence of the partial molar 
thermodynamic properties of infinitely dilute solutions. Estimation of the values of the Krichevskii parameter, Ax,, 
would be helpful in improving the quality of predictions of the infinite dilution partial molar properties of aqueous 
nonelectrolytes at near-critical and supercritical conditions. In this paper, we discuss the methods to obtain values of 
Ax, from various types of experimental information on the thermodynamic properties of binary systems, and use 
these methods to evaluate the Krichevskii parameter for more than 30 aqueous nonelectrolytes. The possibility of 
calculating the values of the Krichevskii parameter for solutes at infinite dilution in water from the van der Waals and 
Peng—Robinson equations of state was also checked. It is shown that the classical cubic equations of state may be 
used to predict the values of Ax, for small nonpolar solutes in water, but grossly overestimate Ax, for larger 
compounds. As an alternative, we propose an empirical correlation between the Krichevskii parameter and the Gibbs 
energy of hydration at 298.15 K and 0.1 MPa for the corresponding compound. This correlation seems to hold for 
solutes with values of the Krichevskii parameter ranging from — 200 to + 200 MPa, and values of the Gibbs energy 
of hydration extending over 100 kJ mol~', and appears to be useful for estimating values of the Krichevskii 
parameter for many solutes of different sizes and polarities, from inert gases to aqueous silica. Predictions, related to 
the sign of the Krichevskii parameter for a number of aqueous inorganic and organic nonelectrolytes, are made, 
which can be verified experimentally. © 2001 Elsevier Science B.V. All rights reserved. 


Keywords: Critical point; Water; Krichevskii parameter; The Gibbs energy of hydration; Aqueous nonelectrolytes; Correlation 











1. Introduction 


There is growing demand for accurate estimates 
of thermodynamic properties of aqueous solutions 
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termination of the partial molar volumes and 
heat capacities of aqueous solutes at near-critical 
conditions [1—3], the high-temperature behavior 
for a vast majority of aqueous solutes has to be 
predicted in the framework of existing semiem- 
pirical models that correlate the thermodynamic 
properties of solutes as functions of state vari- 
ables. A number of such models exist, among 
which the revised Helgeson—Kirkham-—Flowers 
model (the revised HKF) stands out as the 
model with parameters for hundreds of ions and 
nonelectrolytes [4—6]. Recently, several other 
correlating equations were proposed for these 
purposes [7,8]. 

The derivatives of the partial molar Gibbs en- 
ergy, GS, of a solute at infinite dilution diverge 
approaching the critical point of the solvent due 
to the divergence of the compressibility of the 
solvent at these conditions [9—12]. The diverging 
properties include the partial molar enthalpy, 

S, entropy, S$, heat capacity, Cp$, and vol- 
ume, VS. Depending on the nature of solute— 
water interactions, either the ‘+’-type of 
divergence (diverging properties go to infinity), 
or the ‘—’-type of divergence (diverging proper- 
ties go to negative infinity), develops. All elec- 
trolytes in water demonstrate the ‘—’-type of 
divergence, and all aqueous nonpolar nonelec- 














300 400 500 600 700 
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Fig. 1. Partial molar volumes of aqueous methane and boric 
acid at different temperatures and 28 MPa [1,3]. 


trolytes demonstrate the ‘+ ’-type of divergence. 
However, the polar nonelectrolytes may develop 
either ‘+ ’-type, for example, ammonia NH, 
[2,3], or ‘—’-type, as boric acid H;BO, [1]. 

One of the biggest challenges for any model 
for the thermodynamics of aqueous solutions is 
to predict correctly the sign and magnitude of 
the divergence of the partial molar properties of 
infinitely dilute polar nonelectrolytes in the 
vicinity of the critical point of pure water. The 
experimental data, usually available only at 
room temperature, are of little help to deduce 
the near-critical behavior, because of the over- 
whelming contribution of the intrinsic properties 
of the solute to the numerical values of V$ and 
Cp$ at ambient conditions. For instance, meth- 
ane (CH,) and boric acid (H;BO;), two com- 
pounds of about the same size, have very 
similar values of the partial molar volume up to 
500 K (see Fig. 1), and the sign of near-critical 
divergence is obvious only at temperatures 
above 550-600 K. As a rule, there are no ex- 
perimental data for the majority of solutes at 
such high temperatures. If a wrong guess is used 
to deduce the near-critical behavior of a solute, 
then there will be a systematic shift in the pre- 
dicted values of the standard thermodynamic 
functions of this solute at high temperature and 
pressure. For example, if Cp$ for methane is 
used as an approximation of Cp$ for boric acid 
at temperatures above 450-500 K at 28 MPa, 
then the predicted values of the equilibrium con- 
stants for the reaction H;BO,(g)<>H,BO,(aq) 
will be systematically too positive. The differ- 
ence between ‘correct’ and ‘wrong’ values of the 
equilibrium constant In for this reaction are 
plotted in Fig. 2. The magnitude of the error 
increases rapidly with temperature and exceeds 
one In unit in the near-critical region. 

The quality of predictions of the infinite dilu- 
tion high-temperature thermodynamic properties 
of aqueous nonelectrolyte solutions can be im- 
proved given the means to estimate the sign and 
the magnitude of the near-critical divergence of 
partial molar properties for these solutes in wa- 
ter. The search for these means is the main goal 
of this paper. 
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Fig. 2. The error in the In K value for the reaction H,BO,(g)<> 
H,;BO,(aq) due to the wrong approximation of the infinite 
dilution partial molar heat capacity of the solute, see text. 


2. The Krichevskii parameter 


The thermodynamics of infinitely dilute solu- 
tions near the solvent’s critical point are governed 
by the value of the derivative (6P/0X>)).7.y.~o at 
the critical point of the pure solvent [9,11—13], 
which Levelt Sengers [13] has called the 
Krichevskii parameter, A,,: 


- 0X, TVX ,=0 


where P designates pressure, X, represents the 
mole fraction of a solute, V stands for volume, 
and the superscript c indicates that the evaluation 
is done at the critical point of the pure solvent. 
The Krichevskii parameter is finite and has the 
dimension of pressure. For aqueous solutes, the 
Krichevskii parameter is negative for all elec- 
trolytes and for some polar nonelectrolytes in- 
cluding boric acid (H;BO,) and dissolved silica 
(SiO), and positive for most nonelectrolytes. The 
strongly diverging infinite dilution partial molar 
properties of solutes in the vicinity of the critical 
point of water can be represented as the product 
of the Krichevskii parameter and the well known 
thermophysical properties of pure water. Exam- 
ples of these connections are given by [13]: 


5 ~ Ax( ViKr) 
5 ~ Ax( TV }ap) 


Here V? designates the molar volume of pure 
water, x, and &p represent the compressibility 
and the expansion coefficient of pure water. 
These relations are exact at the critical point of 
water and approximate in its vicinity. 

The physical meaning of the Krichevskii 
parameter can be visualized as the measure of the 
change in pressure at the critical point of a sol- 
vent in the imaginary case where a solvent 
molecule is transformed into a solute molecule at 
constant temperature and pressure [14]. In terms 
of statistical mechanics [15,16], the Krichevskii 
parameter for aqueous systems reflects solute— 
water correlations at the critical point of water 
— if the solute—water direct correlations are 
more positive than the water—water direct corre- 
lations, then Ay, <0; if the solute—water direct 
correlations are less than the water—water direct 
correlations, then A,,>0. As the Krichevskii 
parameter determines the sign and the magnitude 
of the near-critical divergence of the partial molar 
properties for solutes, it would be very useful to 
generate estimates of this property for a variety 
of nonelectrolytes to improve the quality of pre- 
dictions of high-temperature thermodynamic 
functions for these solutes in aqueous solutions. 


3. Evaluation of the Krichevskii parameter for 
aqueous nonelectrolytes from experimental data 


There are several ways to evaluate the 
Krichevskii parameter for aqueous nonelec- 
trolytes from published experimental information. 
1. In the first method, the value of the 

Krichevskii parameter is calculated from the 

initial slopes of the critical line of a binary 

mixture [9,12]: 


eh ee -(4)(#) (4) 
. dX3/ cri dT) ,\dX3/ cre 


or the related equation 


4 ol oe -(¥) () 
LAM fees: VAT, FX fee 
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where the subscript CRL stands for the critical 
line, the subscript o indicates the vapor pres- 
sure curve of the pure solvent, and for water 
(dP/dT)S = 0.267 MPa~'. The values of Ax, 
calculated in this way may bear errors up to 
20—25%, owing to experimental uncertainties 
and the often complicated shape of the critical 
line in P—T—X space. 

The second group of methods is based on 
the analysis of the expansion of the Helmholtz 
energy of a mixture at the critical point of the 
solvent [11,12]. Many cited results were first 
derived by Japas and Levelt Sengers [17]. 
When obtained in this way, asymptotic rela- 
tions for properties like Henry’s constants, 
distribution factors, etc. contain a dominating 
term that is the product of the Krichevskii 
parameter and the pure solvent properties. 
These methods are not valid for solutes with 
values of Ay, close to zero, as the product 
term is no longer the dominating term in the 
expansion. This case was first described and 
explained for aqueous D,O by Japas et al. [18]. 
. Values of Ay, can be calculated from the 
near-critical variations of the Henry’s constant 
[17,19], where the following relation must 
hold: 


T n( le A+ B(p — p.) (5) 


where k,, represents Henry’s constant defined 
as ky, = lim of f5/X>) and f, stands for the 
fugacity of a solute; f; and p designate the 
fugacity and liquid-phase density of the pure 
solvent along the gas-liquid saturation curve; 
p, stands for the solvent’s critical density, and 
A and B are constants. The constant B is 
related to the Krichevskii parameter by 


Ax, = Ro:B (6) 


However, it is firmly established [14,20,21] 
that the true asymptotic slope is reached only 
within a few K of the critical point of water, 
where accurate k,, data are virtually absent. 
Therefore, this method is of little practical 
importance. 

. The Krichevskii parameter can be calculated 
from variations in the vapor-—liquid distribu- 


tion factor [17,20,21], K,, which is defined as 
Kp = lim o( Y2/Xo)r5 where Y, and _X, are the 
molar fractions of the solute in the vapor and 
liquid phases, respectively. It was shown 
[17,20] that at near-critical conditions, the fol- 
lowing relation must hold: 

T In Kp =2B(p — p,) (7) 
where B is the same constant appearing in 
Eqs. (5) and (6). Plotting 7 In Ky versus (p — 
p.) allows estimates of B, and, by means of 
Eq. (6), Ax, It was found [20,21] that B is 
nearly constant over the 200 K_ below the 
critical point of pure water. The evaluation is 
facilitated by the fact that at T= 7, the value 
of T In Ky = 0, because Ky = 1. Alvarez et al. 
[21] presented values of B for many gas—water 
systems, which are expected to have uncertain- 
ties of a few per cent. These data can be 
converted to the values of Krichevskii parame- 
ters with uncertainties within 5-10 MPa. This 
is probably the most precise set of Ax, values 
for aqueous solutes. 


. The value of the Krichevskii parameter can be 


evaluated from the initial slope of isothermal 
dew—bubble curves, in P—X coordinates, near 
the solvent critical point [12,22]. As derived by 
Levelt Sengers [12], the expansion of pressure 
on both the dew and bubble sides at the 
critical temperature of the solvent results in 
the expression: 

Px P.+Ay,Xo+°" (8) 
where X, stands for the compositions of either 
‘liquid’ or ‘vapor’ phases, and P. represents 
the critical pressure of a solvent. In other 
words, near the critical point of a solvent, the 
isothermal dew and bubble curves have a joint 
slope equal to the Krichevskii parameter. 
From the practical point of view, to apply this 
method, one should take the average of the 
slopes of dew and bubble curves at tempera- 
tures close to the critical temperature of a 
solvent, as illustrated by Levelt Sengers in [22]. 
For aqueous systems, the coordinates of the 
dew-—bubbles curves are usually presented 
graphically, and estimates of the Krichevskii 
parameter by this method bear considerable 
uncertainties. 
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The third group of methods uses the fact 
that diverging infinite dilution partial molar 
properties of solutes near the critical point of 
a pure solvent can be replaced by the 
product of the Krichevskii parameter and the 
well known properties of the solvent [13]. 


_ Near-critical variations of the partial molar 


volumes V$ provide the means to evaluate 
Ax,. From Eg. (2) one can write 


V3 
i 9 
T> tim (24) si 


The usefulness of this relation follows 
from the empirically established fact that the 
derivative (CP/0X>)y,7.x, =0 for both aqueous 
[23] and nonaqueous [24] systems has practi- 
cally no temperature dependence in the near- 
critical region, at least at reduced 
temperatures (7,=T7/T.) up to 1.05—1.10. 
So, one can plot values of V$/x,V? against 
density and take the value of this function at 
p. as the Krichevskii parameter, as we have 
done in Fig. 3 for benzene. 


. Analogously, the Krichevskii parameter can 


Vo MV1°K;), MPa 


be estimated from the near-critical variations 
of the partial molar enthalpy, H$. From Eq. 
(3), one can write 


3 
A= 10 
‘ged lim, (a5) ( ) 





Critical density 











250 300 350 
Density, kg/m?* 


‘“‘However, as the partial molar enthalpy is a 
relative property, containing an arbitrary ref- 
erence constant, this relation cannot be used 
for quantitative determination Of Ax, ” [25] 

. In principle, the near-critical variations of 
the infinite dilution partial molar heat capac- 
ities, Cp5, can be used to evaluate the values 
of the Krichevskii parameter by means of the 
relation Cp$=(@H3/0T)p. However, our at- 
tempt to treat Cp$ data from Hnedkovsky 
and Wood [3] in this way showed that the 
resulting function exhibited considerable scat- 
tering, presumably due to very strong tem- 
perature variations in the term (0%p/0T)>, at 
near-critical conditions, where the uncertainty 
of up to 0.7 K in the temperature control of 
the experiments by Hnedkovsky and Wood 
[3] might be too large to allow quantitative 
determination of the Krichevskii parameter 
from Cp§ data. 

We used these relations to evaluate values of 
the Krichevskii parameter for a large number of 
aqueous nonelectrolytes. The calculation of Ax, 
from the initial slopes of the critical line re- 
quires experimental results at very low solute 
concentrations, particularly for binaries with 
complicated critical curve shapes. So, the coordi- 
nates of the critical line of a binary mixture 
were used only if there were data at the critical 
mole fraction, X,, less than 0.05 (for the system 
ethanol—water, the minimal critical mole frac- 
tion is about 0.07, and it was used because this 
system had a monotonous (without extremes) 
critical line [27]). Typical uncertainties in values 
of Ax, are expected to be within + 25-40 MPa 
if evaluated from the initial slopes of the critical 
lines or the dew—bubble curves, and +10 MPa 
from V3 data. The necessary thermophysical 
properties of pure water in the near-critical re- 
gion were calculated from the Hill equation of 
state [28]. The resulting values of the 
Krichevskii parameters for 32 aqueous nonelec- 
trolytes are given in Table 1, together with val- 
ues of the Gibbs energy of hydration, A,,G°, in 
kJ mol~', at 298.15 K, 0.1 MPa (corresponding 


Fig. 3. Illustration of the use of Eq. (9) in the text to evaluate 
the value of the Krichevskii parameter for benzene in water 


from V$ using data from [26]. 


to the process of transferring 1 mole of a solute 
from an ideal gas state to an ideal one molal 
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Table | 


The evaluated values of the Krichevskii parameter for aqueous nonelectrolytes 








Nonelectrolyte 


A,G°, 
(kJ mol~') 


CRL 


DISTR* 


Dew-—Bubble From V$ 


From H$ 








n-Dodecane, C,H. 
Tetrafluoromethane, CF, 
Helium, He 
Neon, Ne 
n-Heptane, C,H, 
Nitrogen, N, 
n-Hexane, C,H,,4 
Hydrogen, H, 
n-Pentane, C;H,> 
Oxygen, O, 
Argon, Ar 
Methane, CH, 
Propane, C,H, 
Ethane, C,H, 
Krypton, Kr 
Xenon, Xe 
Ethene, C,H, 
Cyclohexane, C,H,> 
Propene, C3H, 
Carbon dioxide, CO, 
Hydrogen sulfide, HS 
Toluene, C,H, 
Benzene, C,H, 
Ammonia, NH; 
Cyclohexylamine, C,H,3;N 
Ethanol, C,H,O 
Piperidine, C;H,,N 
Quinuclidine, C;H,;N 
Water, H,O 
Heavy water, D,O 
2-Amino-2-methyl-1-propanol, 
C,H,,NO 
Boric acid, H;BO, 
Silica, H,SiO, 








2 ar° 
20.979 
19.45¢ 
19.08° 
18.85° 
18.20° 
18.12° 
i. 
17.40° 
16.55¢ 
16.29¢ 
16.294 
16.09° 
15.57° 
14.849 
13.454 
13.25¢ 
12.99° 
12.87° 
8.38¢ 
5.66P 
4.56” 
4.29» 
— 10.05" 
— 12.48 
— 13.00° 
— 13.48 
— 14.5° 
— 18.55" 
— 18.90" 
— 26.3° 


— 36.9" 
— 69.3" 


Lt ar 
210 + 20° 


234+ 40° 
180® 
151+15° 
150° 
181 + 30° 
140® 
130° 
150° 
160 + 20! 


1s2 + 29° 
132 4%25° 
1008 


66 + 254 


$2427 


(0) 
0.48” 


40 
57 + 20° 


sit oer 
49 + 20° 
(0) 
0.47° 
—6 + 20° 


— 67.5” 


167 + 40° 


197 + 40" 170 + 10! 


180 + 40° 


188 + 15) 
163 + 10* 
153 + 40! 
194 + 40™ 


166 + 15) 
176 + 15) 
130 + 15! 


130 + 10* 
105 + 10 
105 + 10! 
95 + 10! 
52 + 10 


— 67 + 10% 
— 184 + 20” 


190 + 15! 


160 + 15! 


125 + 15! 
115 + 15! 





Eq. (17) 


181 
175 
169 
168 
167 
164 
164 
162 
161 
157 
156 
164 
155 
153 
150 
146 
144 
142 
142 
123 
112 
108 
107 
48 
38 
35 
34 
29 
12 
1] 
—19 


— 63 
— 196 


* Alvarez et al. [21], unless otherwise indicated; >Plyasunov and Shock [47]; “Stevenson et al. [48];“Wilhelm et al. [49]: “Smits et 
al. [50]; "de Loos et al. [51]; As given in Harvey et al. [20]; "de Loos et al. [52]; ‘Degrange [26]; /Biggerstaff and Wood [53]; 
‘Hnedkovsky et al. [2]; 'de Loos et al. [54]; ™Danneil et al. [55]; "Brollos et al. [56] as given in Hicks and Young [57]; “Sanchez and 
Lentz [58]; PCarroll and Mather [59]; 4Rebert and Kay [60]; "Kawazuishi and Prausnitz [61]; ‘Lewis and Wetton [62]; ‘Griswold et 
al. [27]; “See text; ‘Japas et al. [18]; “Kukuljan et al. [63]; “Hnedkovsky et al. [1]; ‘Fournier and Potter [32]. 


aqueous solution). The third, fourth, fifth, and 
sixth columns indicate the type of experimental 
data used to evaluate Ax, — the initial slopes of 


the critical line (CRL), the vapor—liquid distribu- 
tion coefficients (DISTR), the isothermal dew— 


bubble curves (DEW-B), the infinite dilution 
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partial molar volumes (V$), respectively. The last 
column gives the values of the Krichevskii 
parameter calculated using the correlating equa- 
tion (Eq. (17)), see below. 

The Gibbs energy of hydration of most solutes 
at 298.15 K, 0.1 MPa are taken from literature 
compilations. The standard state adopted for an 
aqueous nonelectrolyte is the hypothetical state of 
the pure substance in water of unit activity at 
molality equal to one at any temperature and 
pressure; the standard state adopted for a gas is 
the ideal gas of unit fugacity at any temperature 
and pressure of 0.1 MPa. The Gibbs energy of 
hydration of pure water is equal to the Gibbs 
energy of vaporization of liquid water, A,,,G°, 
with the opposite sign, recalculated into the mo- 
lality concentration scale, A,G°=A,,,G° — 
RT 1n(1000/M,,) = — 8.59 —9.96 = — 18.55 kJ 
mol~', where M,, = 18.0152, and g mol! is the 
molar mass of water. The Gibbs energy of vapor- 
ization of liquid water is taken from the CO- 
DATA recommendations [29]. The Gibbs energy 
of hydration of boric acid is equal to the sum of 
the standard Gibbs energy of solution of crys- 
talline boric acid, A,,,G°, and the Gibbs energy of 
sublimation of crystalline boric acid, A,,,G°, 
taken with the opposite sign, A,G° = A,,,G° — 
AgwiG° = (0.392 + 0.020) -- (37.3 + 2.3) = — 36.9 
+2.3 kJ mol~!. The numerical values of the 
Gibbs energy of solution and the Gibbs energy of 
sublimation of crystalline boric acid are taken 
from [29,30], respectively. The Gibbs energy of 
hydration of D,O in the molality concentration 
scale was calculated by means of A,G° = RT 
In(¥P./P®) + RT In y** — RT 1n(1000/M,,) = 
— 18.90 kJ mol~', where P,=0.00271 MPa 
stands for the saturated vapor pressure of heavy 
water [31], ¥% ~ 1.0 represents the fugacity coeffi- 
cient of pure D,O; y** = 1.0 is the activity coeffi- 
cient of D,O in its (nearly ideal, see [18]) binary 
mixture with H,O. For dissolved silica, SiO;, the 
infinite dilution partial molar volumes may be 
evaluated from the equation that correlates the 
solubility of quartz in water over a wide range of 
temperatures and densities [32]. Note that the 
value of the Krichevskii parameter does not de- 
pend on the extent of hydration of silica in 


aqueous solution. However, the Gibbs energies of 
hydration of various hydrates of SiO, are very 
different, —69.3+14 kJ mol~' for H,SiO,; — 
315.0 + 26 kJ mol~! for H,SiO;; — 548.9+9 kJ 
mol~' for SiO,, all at 298.15 K and 0.1 MPa. 
These values are calculated by combining the 
Gibbs energy of formation of hydrates of SiO, in 
the gaseous state as obtained by ab initio calcula- 
tions in [33], the partial molar Gibbs energy of 
formation of SiO, in aqueous solution as evalu- 
ated in [5], and the Gibbs energy of formation of 
H,O from the CODATA recommendations [29]. 
As it is widely assumed that the predominant 
form of silica in aqueous solution is H,SiO, [34], 
we tabulated in Table 1 the Gibbs energy of 
hydration of this species. 


4. Correlations to estimate the Krichevskii 
parameter for aqueous nonelectrolytes 


4.1. Calculations of the Krichevskii parameter in 
the framework of the van der Waals and the 
Peng— Robinson equations of state 


We examined the possibility of calculating the 
Krichevskii parameter, A,,, within the framework 
of the classical equations of state (EoS) for wa- 
ter—nonelectrolyte mixtures. This can be done 
through the link between A,, and the infinite 
dilution solute—solvent direct correlation function 
integral at the critical point of a pure solvent, 
Ci*(T.=1, p,=1), where values of both the 
reduced temperature, 7., and the reduced density, 
P,, are equal to one. O’Connell [35,36] obtained 
the general relation between the infinite dilution 
partial molar volume of a solute and C%*: 


where C¥,* stands for the infinite dilution solute— 
solvent direct correlation function integral at arbi- 
trary temperature and density. Comparing Eqs. 
(11) and (2) one finds the relation between Ax, 
and Cy,*(T, = 1, p, = 1): 


= ( a Ci,*(T, = I, Pr= 1))RT, 
Vi(T.. Po) 


(12) 





Ax, 
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Petsche and Debenedetti [37] analyzed the effect 
of the difference in size, energy and shape between 
solute and solvent on the sign of the near-critical 
divergence of the partial molar properties. For the 
van der Waals EoS these authors derived the 
following relation to calculate C¥#,*(T.=1, p,= 


1): 
Ci2*(T, = 1, p, = 1) 


_3Q +7) , 90°51 ~ ki) 
4 4 


where a=a,/a,, B =b,/b,, a,, b,, a5, b, are the 
parameters of the van der Waals EoS for solvent 
(water) and solute, respectively; and k,, is a mix- 
ture-specific parameter appearing in the mixing 
rule for the a parameter of a mixture (see Ap- 
pendix A): 


> = (a,az)(1 — ky») (14) 


We used Eqs. (12)—(14) to predict the value of 
the Krichevskii parameter for some solutes given 
in Table 1. The necessary values of k,, were 
evaluated from the second cross virial coefficient 
between water and solute, B,,, using the Hayden— 
O’Connell correlation [38] at one temperature ar- 
bitrarily chosen equal to the critical temperature 
of pure water, by means of ([39], p. 33): 


=] (13) 





ky.=1- (0.5(b, + by — B,»)| (15) 


ae 
(a 14>)? 


Besides applying the procedure described by 
Petsche and Debenedetti [37], we derived the rela- 
tion to calculate C*#,*(T.=1, p,=1), from the 
Peng—Robinson [40,41] EoS, which is more prac- 
tical for engineering applications (see Appendix 
A): 

Ci.*(7, = 1, p, = 1) 
4.578) 
(1 — 2) 
ae 22,271 — 2) 
(1 + 22, — 25)” 
where 22, = 1.489 and Q, ~ 0.2534. Again Eq. (15) 
was used to evaluate k,, values. 


Except the case of D,O using the van der Waals 
EoS, both equations correctly predict the positive 


22,0°°(1 ps ky>) 
1 +20, — Q? 





(16) 





sign of the near-critical divergence for solutes 
tested (comparisons for some polar compounds 
cannot be made due to lack of the necessary 
critical parameters). For small solutes (inert and 
other nonpolar gases, CO,, NH;) the agreement 
with experimental results is often quite satisfac- 
tory, within +20 MPa. However, for larger so- 
lutes (benzene, n-hexane), the calculated values 
are grossly overestimated; the deviations being 
100 MPa and more. The reasons may be an 
inadequacy of the simple mixing rules for compo- 
nents of very different sizes or a general inappli- 
cability of cubic equations of state to water, a 
substance with a very low critical compressibility 
factor. Another complicating circumstance is the 
strong sensitivity of results to the value of the 
mixture-specific parameter k,5, particularly for 
polar compounds. As a result, we expect that the 
cubic equations of state may be used only for 
estimating values of the Krichevskii parameter for 
mixtures of water with small nonpolar solutes. 


4.2. An empirical correlation between the 
Krichevskii parameter and the Gibbs energy of a 
solute at 298.15 K, 0.1 MPa 


Hnedkovsky [42] noticed the correlation be- 
tween the VS values of CH,, CO,, HS, and NH, 
at extrema of isothermal compressibilities at near- 
critical conditions and the Gibbs energy of hydra- 
tion at 298.15 K, 0.1 MPa of these compounds. 
As V3 values in the near-critical region are pro- 
portional to the Krichevskii parameter, see Eq. 
(2), this finding can be reformulated as a correla- 
tion between Ax, and A,G° at 298.15 K, 0.1 MPa. 

In Fig. 4, values of the Krichevskii parameter 
for solutes are plotted versus the Gibbs energy of 
hydration, all from Table 1. The circles, triangles 
and squares represent Ax, values evaluated from 
the initial slopes of the critical lines and from the 
isothermal dew—bubble curves, from vapor-liq- 
uid distribution constants, and from VS data, 
respectively. The size of symbols is close to the 
expected uncertainties in the values of the 
Krichevskii parameter. The diamond gives the 
value for water as a solute in water as a solvent. 
The horizontal line in Fig. 4 shows the value of 
Ax, =0, the vertical line locates the value of the 
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Gibbs energy of hydration of H,O. The linear 
correlation in this plot is given by: 


although the reasons for its success remain ob- 
scure. The calculated value of Ax, is in MPa, the 
Gibbs energy of hydration is in kJ mol~'. This 
linear correlation seems to be valid for com- 
pounds with values of the Gibbs energy of hydra- 
tion between — 80 and 25 kJ mol! at 298.15 K, 
0.1 MPa, with values of the ‘excluded’ volumes 
(taken as the molar volume at the normal boiling 
point of the compound) ranging from 19 (water) 
to at least 150 cm? mol~! (n-hexane). We expect 
that the uncertainty is about 20-25 MPa in the 
predicted values of the Krichevskii parameter. Eq. 
(17) does not give an exact result for water, and 
attempts to force it to go through the value of 
Ax, =0 for water degrade the performance of the 
correlation for solutes with large absolute values 
of the Gibbs energy of hydration. Values of the 
Krichevskii parameter calculated by means of Eq. 
(17) are given in the last column of Table 1. 


Several points should be noticed. First, this 
correlation is empirical. Indeed, although both the 
Krichevskii parameter and the Gibbs energy of 














-40 20 0 20 40 
AnG°(298 K, 0.1 MPa), ku/mol 


Fig. 4. The correlation for aqueous nonelectrolytes between 
the values of the Krichevskii parameter and the Gibbs energy 
of hydration at 298.15 K, 0.1 MPa. The line represents Eq. 
(17). Different symbols represent different sources of experi- 
mental data used to evaluate A,,, see text. The size of these 
symbols is close to the uncertainty of A,,. The diamond gives 
the value for water. 


hydration at 298.15 K, 0.1 MPa reflect the 
strength of solute—water interactions, they refer 
to different properties at different temperatures 
and densities. For example, the relative contribu- 
tions of cavity formation work and solute—water 
interactions are expected to be quite different at 
room temperature and high density than at criti- 
cal values of the temperature and density of pure 
water. Second, Eq. (17) predicts the change of the 
sign of the Krichevskii parameter at a single value 
of the Gibbs energy of hydration, close to the 
value of the Gibbs energy of hydration of water 
as a solute. However, as solutes differ in polari- 
ties, shapes, number of groups able to dissociate, 
etc., it is expected that no single empirical correla- 
tion can universally locate the position of the sign 
change for Ax,,. So, we call for caution when 
using Eq. (17) for compounds with Gibbs energies 
of hydration that are within + 3-5 kJ mol~! of 
the corresponding value for water. Third, scaled 
particle theory [43] shows that there is a strong 
dependence of the Gibbs energy of hydration at 
high densities on the size of a solute. As a conse- 
quence, the proposed relation may not be valid 
for solutes of very large size (much larger than 
n-hexane), although this remains to be tested. 

A number of conclusions, which can be verified 
experimentally, can be drawn using literature 
compilations of the Gibbs energy of hydration of 
inorganic [44] and organic compounds [45]. Most 
inorganic ‘gases’ are expected to have positive 
values of the Krichevskii parameter (some exam- 
ples are SO,, H,Se, O;, Cl, Br,, I, IC, IBr, 
HCIO, HN;, HNO,, PH,, NOCI, and Hg), a few 
are expected to have negative values of Ax,, 
H,0,, N,H,, SiF,, and neutral ionic pairs like 
HCl and HNO,. Among the numerous organic 
compounds most hydrocarbons are likely to have 
positive values of the Krichevskii parameter with 
the possible exceptions of polyaromatic com- 
pounds with more than four condensed rings. 
Monofunctional, nonaromatic compounds includ- 
ing alcohols, ketones, amines, ethers, esters, nitro- 
compounds, and nitriles are expected to have 
positive values of the Krichevskii parameter. 
Some organic compounds, such as phenol and 
monocarboxylic acids, have values of A,G° at 
298.15 K, 0.1 MPa that are close to that of H,O. 
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It is expected that compounds derived by addition 
of nonpolar groups like CH, and CH, (which 
increase the Gibbs energy of hydration) to these 
‘border’ solutes, will have more positive values of 
Ax, compared with ‘parent’ compounds. Like- 
wise, the addition of polar groups, such as OH, 
CO, COOH, NH>, CONH,, or benzene ring, will 
decrease values of the Krichevskii parameter. So, 
we expect that polycarboxylic acids, hydroxy 
acids, amino acids, and benzenecarboxylic acids 
will all have Ax, <0. Recent experimental deter- 
minations of apparent molar volumes of the 
aqueous amino acids a-alanine, B-alanine and 
proline at elevated temperatures by Clarke and 
Tremaine [46] showed that values of V$ pass 
through maxima at 400—450 K and then decrease, 
confirming the negative values of the Krichevskii 
parameter for these amino acids. Amides, polyols 
and carbohydrates have more negative values of 
the Gibbs energy of hydration than H,O, and we 
expect that these compounds also have Ax, <0. 
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Appendix A. The statement for C7,* for the 
Peng—Robinson EoS 


We followed the same procedure that was used 
by Petsche and Debenedetti [37] for the van der 
Waals EoS. The Peng-—Robinson EoS is written 
as 
2 


pweT ap 

P= ae > 5 
l—pb 1+2bp —b-p? 
where p=N/V and N=N,+N,. The following 
mixing rules are used, a= Xia, +2X,Xa,.+ 
X3a,, b= X,b,+ Xob2, where a,> = (a,a>)°°(1 — 


(Al) 





k,,) and X,=N,/N. Index 1 refers to a solvent, 
index 2 refers to a solute. The expression for 6 is 
given by: 


‘ oP 
o= lim x(: 
x20 ON, TLV.N | 


_ pKT{l + p(bz— 5,)) | 2d) 
(i ~ pir 1 + 2b,p — bip* 
2a,b,p[l — b2p] 
(1 + 2b,p — bip*) 


Finally, we obtain for the infinite dilution direct 
correlation function integral 


| + p(b, —b,) 
(l= pb, 
‘ 2a1> 
pKT(1 + 2b,p — bip’) 
2a,b,p[1 — bp] 


~ pxT(1 + 2b,p — b?p?P° 
(A3) 








(A2) 











The next step is to express C*,* in terms of 
reduced properties of the pure solvent. For the 
Peng—Robinson EoS a= Q,(kT./p f(T), 6 = Q,/ 
P., Where Q, x 1.489 and Q, ~ 0.2534. The tem- 
perature dependent function is given’ by 
f(T) =[1+m(1 — T°°))?, and the parameter m is 
m = 0.37464 + 1.54226 — 0.26992m", where @ 
stands for the Pitzer acentric factor of a com- 
pound. We define y = b,/b,, « =a,/a,, and obtain 





22,p,0°°(1 rn k1>) 
T,(1 + 22,9, — Q¢p7) 
= 22,2479 71 — Qyp,) 
TAI + 225: tid Qpi) 





(A4) 





At the critical point of a pure solvent we have 


Ci,*(T= 1, p, = 1) 
tae) 
(1 — Q,)° 
22,247 (1 — 2.) 
~ (1+ 2, — 227 


22,,2°°(1 — ky2) 
1+ 22, — 23 





(AS) 
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Abstract 


A predictive model based on the Krichevsky—Ilinskaya equation (KIE) combined with Peng Robinson EOS and the 
modified Huron—Vidal mixing rules has been proposed to estimate hydrocarbon solubility in both the vapour and 
liquid phases. For solubility in the liquid phase, the proposed model estimates the partial molar volume from an 
equation of state with Huron—Vidal mixing rules whereby a modified UNIFAC model is used to predict the 
asymmetric activity coefficient. The hydrocarbon data considered for modeling were non-polar liquid hydrocarbons 
such as n-alkanes and aromatic benzene compounds. For the hydrocarbons studied, the proposed model satisfactorily 
predicted the supercritical carbon dioxide solubility in the vapour phase. It was comparable to that correlated using 
equations of state with fitted interaction parameter applied to both phases. However, the liquid phase were less 
satisfactorily predicted with an AAD ranging from 5 to 24%. Partial molar volume at infinite dilution was found to 
be a critical input parameter as changing this parameter improved the error deviation between 4.1 and 10%. © 2001 


Elsevier Science B.V. All rights reserved. 
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Nomenclature 


Parameter in EOS [m°atm/kgmo!?] 
Interaction parameter in UNIFAC 
Interaction parameter in UNIFAC 
Interaction parameter in UNIFAC 
Parameter in EOS [m*/kgmol] 
Constant in Eq. (9) 
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Fugacity [atm] 


Excess Gibbs free energy [atm m°*/kgmol] 
Henry’s constant of a solute gas in a solvent [atm] 


Number of moles or carbon number 


Pressure [atm] 


Universal gas constant [m*atm/kgmol K] 


Molecular volume parameter in UNIFAC 


Temperature [K] 


Partial molar volume at infinite dilution [m*/kgmol] 


Liquid mole fraction 
Gas mole fraction 


Greek letters 
y Activity coefficient 


y Unsymmetric activity coefficient 
wy) Fugacity coefficient 


Subscripts 
i Component (i ) 
m Mixture 


Superscripts 
G Gas phase 
L Liquid phase 











1. Introduction 


The study of supercritical carbon dioxide 
(SCCO,)-hydrocarbon systems is necessary be- 
cause of its increasing application in enhanced oil 
recovery by miscible gas displacement and as an 
extraction solvent. There is a lot of solubility data 
available for a variety of solute-solvent systems. 
However, modeling of such systems, especially for 
larger asymmetrical hydrocarbons is difficult. A 
convenient approach to model solubilities at su- 
percritical conditions is the use of equations of 
State applied to both phases with various mixing 
rules. Equations of state have been used with van 
der Waals mixing rules (VaWMR) to model sys- 
tems of SCCO, and crude oil [1], alkenes and 
m-xylene [2], and aromatic hydrocarbons [3]. 

Much work has been done in modeling high 
pressure vapor—liquid and gas—liquid phase equi- 
libria and many prediction methods have been 
proposed [4—8]. Since the proposition of Vidal [4] 


to combine an excess Gibbs free energy (G*) 
model with an equation of state (EOS), where the 
G® from an EOS evaluated at infinite pressure 
was equated with that of an existing G* model, 
the method has been receiving much attention and 
attempts have been made to improve its predictive 
accuracy. Huron and Vidal [5] improved the 
method to allow simple EOSs to describe complex 
systems. Michelsen [6] proposed the modified 
Huron-— Vidal (MHV) mixing rule by deriving GF 
from an EOS at zero pressure; while Boukouvalas 
and co-workers [7] proposed a model named 
LCVM, which is a linear combination of the 
Vidal and MHV mixing rule. Zhong and Ma- 
suoka [8] indicated that the predictive accuracy of 
MHV mixing rule when applied to highly asym- 
metric systems such as gas/large alkane systems, is 
not satisfactory and becomes poor with the in- 
crease in the alkane carbon number. They also 
noted that although the LCVM model has good 
predictive accuracy for asymmetric systems, it 
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lacks theoretical justification. They proposed a 
new mixing rule, MRI, by introducing a correc- 
tion factor, f, into the MHV mixing rule to im- 
prove the predictive accuracy for highly 
asymmetric gas-large alkane systems. The factor, 
f, is almost a linear function of carbon number of 
heavy alkanes (n>10) and equals to zero for 
n< 10. They tested the method for hydrocarbons/ 
CO,/C,H./CH, systems at pressures up to 10 
MPa and temperatures up to 373 K and for 
alkanes of carbon number up to 36. Similar inves- 
tigations were performed for mixtures of CO, H.,, 
C,H, with large alkanes [9] and for CO, and 
CyoHy2, Cr2H 46, Cr4H5q and CrgHsg alkanes for 
pressures up to 41 MPa [10]. However, MRI 
mixing rule is a correlation that is valid only for 
gas/large alkanes and the factor f is an optimized 
parameter that depends on the type of gas/alkane 
system. 

As the solubility of supercritical fluids in liquids 
is a special case of vapor—liquid equilibria, where 
the gas cannot exist as a liquid at the conditions 
of the solution, improved predictive methods 
based on Henry’s law may be desirable. Equations 
such as that of Krichevsky—Ilinskay (KIE) for the 
prediction of gas solubilities at high pressures has 
been used [11]. For low to moderate gas solubili- 
ties, the predictability of KIE is reasonably good 
[11,12]. The solubility of SCCO, in hydrocarbons 
is very high (could be more than 90 mol% at 
pressures up to 35 MPa). For solution tempera- 
tures not far removed from the critical tempera- 
ture of these solvents, the solution is sufficiently 
expanded and compressible and it is unlikely that 
the partial molar volume becomes independent of 
pressure and composition. 

In this study a predictive model for the solubil- 
ity of SCCO, in hydrocarbons free of any fitted 
parameters has been proposed. The fugacity of 
SCCO, in the liquid phase was calculated by 
estimating the partial molar volume at infinite 
dilution (63°) for the KIE equation due to the 
compressibility of the SCCO,/hydrocarbon solu- 
tion. In addition, the modified UNIFAC model 
[13] was used to predict the unsymmetric activity 
coefficient. The model utilizes the Huron and 
Vidal [5] approach to predict the energy parame- 
ters of the EOS. The predicted values obtained by 


the proposed model were compared with that 
obtained using the Peng Robinson (PR) EOS with 
the vaWMR applied to both phases. For estima- 
tion of binary mixture data using the PR EOS, 
knowledge of interaction parameters was neces- 
sary. This was obtained from literature. 


1.1. Thermodynamic model 


For binary mixtures, the KIE 1s: 


+f 22 ima (1 
n y3, 
Py RT ° iy 2 ) 


where subscript 1 stands for hydrocarbon sol- 
vent and 2 for SCCO, solute. The variation of 55 
with pressure has been included in the KIE ex- 
pression and is given by: 
ee ee 
n—=In Ay) + Inyo + 


x5 


bz (P — Pi) 


RT (2) 





The Henry’s law constant in Eq. (1) and Eq. (2) 
is evaluated at the saturation pressure of the 
solvent from an EOS by [11]: 


HY) = PTS. (3) 


The combinatorial part of the activity coeffi- 
cient,y%, of SCCO, normalized according to the 
unsymmetric convention is determined from the 
modified UNIFAC model and is given by [13]: 


In 9° = In(r79)/(x5r 3? + x17?) 


— (73/3) /(xorZ? + xyrq?) + 37/7”. (4) 


The temperature dependent-interaction parame- 
ters of UNIFAC are given by [13]: 


Aj = ji) + Aji 2(T — To) + 43 TInT/T) + T— To}, 
(5) 


where 7, is an arbitrary reference temperature, 


taken as 298.15 K. 
From an EOS, the partial molar volume at 
infinite dilution is given from the exact relation: 


jms (Gelondese “= 
(OP/o V) T.nj.n; n,=0 
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As for component 1, the fugacity f” is evaluated 
from the following relation: 


b(P — Pi) 
i: eee 


The fugacity of both components in the gas phase 
is expressed as: 


fi = bP y;P. (8) 
The energy parameter expression for the PR EOS 
is given by: 


\ x;a;B 
RT) ~ 25 aRTD, - 
2,/2RTb /2RTb, 


m 


(7) 





Ff = yx, P" exp 





where B is a constant pets as: 
(22~,/2) 


Q+/5 


The volume parameter of the equation of state is 
given by vdW type mixing rules as: 


B= 


bm => bX; (10) 


In this work the excess Gibbs free energy was 
evaluated from the modified UNIFAC method 
[13]. 


2. Results and discussion 


The solubility SCCO, in selected hydrocarbons 
(obtained from literature) in both the gas (y,) and 
liquid (x,) phases was modeled using the PR 
EOS + vdW mixing rule, the PR EOS + HV mix- 
ing rule and the proposed KIE+ PR EOS + HV 
mixing rule methods. The average absolute devia- 
tion (AAD) obtained for the different methods 
and hydrocarbon data reference are shown in 
Table 1. It can be seen that all three methods 
generally give good results for the gas phase solu- 
bility. This is not surprising, especially for the PR 
EOS + vdWMR due to the inclusion of adjustable 
parameters. The PR EOS+ HV estimation did 
not converge for the larger branched alkanes. In 
order to avoid the use of interaction parameters 
and the problem of convergence, the proposed 
KIE + PR EOS + HV model may be a desirable 
choice for predicting the gas phase solubility data 


where a grand average AAD of 1.7% was 
obtained. 

For estimation of CO, solubility in the liquid 
phase, the PR EOS + vdaWMR method was found 
to be superior compared to the other two meth- 
ods, however it is again subjected to prior knowl- 
edge of interaction parameters. The PR 
EOS + HV is unsatisfactory due to deviations and 
non-convergence for the larger branched alkanes. 
The KIE+PR  EOS+HV_ method (with 
o¥calculated using PR EOS) on the other hand 
produced relatively large deviation (grand average 
AAD of 11.4%) and _ non-convergence for 
branched alkanes. This is expected to be primarily 
due to unsatisfactory value obtained for 03 from 
Eq. (6) which was determined by the PR EOS. 
Studies [23,24] have shown that an EOS such as 
PR EOS over predicts, and is always semi quanti- 
tative in the prediction of 05 of gases in liquids. 
There is rarely an EOS that is sufficiently accurate 
for calculating the partial molar volume at infinite 
dilution [11]. Liu and Macedo [25] in their study 
on models for infinite dilution partial molar vol- 
ume of solutes in supercritical fluids concluded 
that no available model is recommended for 
quantitative predictions of infinite dilution partial 
molar volume of solutes in supercritical fluids. 

In order to study the effect of 03 (and hence 
term 2 on the RHS of Eq. (2)) on the errors 
obtained in the liquid phase, the value of 6 was 
arbitrarily varied (assuming independent of pres- 
sure). This is demonstrated in Fig. 1, where the 
AAD of CO, solubility in the liquid phase is 
plotted against arbitrary values for oy for the 
CO, +n-hexane system at 353 K. It can be seen 
that the AAD FOR x, decreases as 05 increases, 
reaches a minimum and then increases with in- 
crease in 63. This clearly shows that large errors 
can result over a small change in 53 and accurate 
determination of this property is important. 

The effect of this parameter for the prediction 
of SCCO, solubility in the liquid phase was evalu- 
ated by using values of 0;°calculated by PR EOS 
and that optimised into the proposed KIE+ PR 
EOS + HV model. This is shown in Fig. 2 for the 
CO,—n-hexane system at 353 K. For this system 
the calculated value of 03 using Eq. (6) with PR 
EOS is 0.071 m?/kgmol. As this value is decreased 
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to an optimum value of 0.03 m?/kgmol the pre- 2.5%. Similar behavior was also noted for most 
dicted solubility approaches the experimental other systems studied. One example the solubility 
value and the AAD of x, reduced from 15 to plot for the CO,—n-heptane system at 352 K 





AAD in x, (%) 


— best fit line 











0.06 


v, (m*/kgmol) 


Fig. 1. Deviation (%) of the liquid phase solubility with varying 5 of CO, for the CO, + n-hexane system at 353 K. 





0.85 


CO, solubility in x, (mol frac.) 


— best fit lines 





T , © t T , FF 7 T 7 T 








P (atm) 


Fig. 2. Solubility of CO, in n-hexane at 353 K: (@) Gupta [14]; (O) 6% =0.071, estimated using PR EOS; (¥) 6 =0.05; (V) 
55 = 0.03. 
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CO, solubility in x, (mol frac.) 


— best fit lines 











P (atm) 


Fig. 3. Solubility of CO, in n-heptane at 352 K: (@) Li et al.; (O) 67 = 0.0663, estimated using PR EOS; (W) 63% =0.04; (V) 
55° = 0.0365. 


Table | 
Deviation of CO, solubility in the vapour and liquid phase obtained using different approaches 











Solvent Data points AAD (%)? Reference 








PR EOS+vdWMR_ PR EOS+HV KIE+PR EOS+HV 











y2 








Hexane 1.3 
Heptane 1.6 
Octane 0.3 
Decane 0.6 
Toluene 10 0.6 
m-Xylene 14 0.8 
o-Xylene 17 0.3 
p-Xylene 18 0.4 
Benzene 15 2.0 
Tetradecane 12 1.3 
Hexadecane 61 0.9 











®@ AAD% = (100/NP) E |(Xexp—Xcat)/Xexpl: 
> Used 5% calculated by PR EOS. 

© Used optimised 05°. 

4 No convergence. 


which is shown in Fig. 3. The solubility data for solubility data is better predicted using the opti- 
x, predicted using the KIE + PR EOS + HV with mised 0% with a grand average AAD of 6.4% 
optimised value of 6% is shown in Table 1. The (excluding branched alkanes). Large improvement 
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AAD of x, decreased from 20.4% at a calculated 
0% of 0.052 m*/kgmol to 4.5% at an optimised o¥ 
value of 0.001 m*/kgmol. The calculated and opti- 
mum values of 03 used in the proposed model are 
shown in Table 2. Generally, reducing the value 
of the calculated 63 improved the prediction and 
reproduction of the literature experimental data 
from 13 to 98%. It has to be acknowledged that 
the calculated 63 values show a large deviation 
from the optimised calculated 65 values and is 
characteristic of an adjustable parameter. 

It can be seen from Table | that the estimation 
method proposed (using optimised value of 03°) is 
able to satisfactorily correlate the x, solubility 
data of SCCO, in most hydrocarbons studied. In 
addition, no interaction parameter was required 
and all the data converged. However, estimation 
of the larger branched alkanes is poor. This was 
also obtained by other investigators [8] using dif- 
ferent models. They proposed new mixing rules 
incorporating a correction factor into’ the 
Michelsen mixing rules [26] to account for large 
alkanes. It is evident from this analysis that the 
estimation of the solubility of SCCO, in liquid 
hydrocarbons will depend on the accuracy associ- 
ated with the prediction of 05. 

The calculation of the Henry’s law constant at 
infinite dilution from equations of state is well 
established. Prausnitz et al. [11] presents various 
examples in which the calculated value of Henry’s 
law constant from an EOS matches very well the 


Table 2 
Calculated and optimised values of 0+ 








Compound Temperature (K) 


Calculated by PR EOS 








0.0710 
0.0663 
0.0623 
0.0590 
0.0530 
0.0624 
0.0480 
0.0498 
0.0520 
0.0792 
0.0714 


Hexane 
Heptane 
Octane 
Decane 
Toluene 
m-Xylene 
o-Xylene 
p-Xylene 
Benzene 
Tetradecane 
Hexadecane 





ox (m?/kgmol) 


experimental data for all systems considered. 
Thus the calculation of Henry’s law constant from 
an EOS using Eq. (3) is fairly accurate and the 
prediction using Eq. (2) is expected to be insensi- 
tive to this parameter. Moreover, the accuracy 
associated with the calculation of y# in Eq. (2) is 
dependant on the calculation of the activity coeffi- 
cients from the modified UNIFAC [13]. 


3. Conclusions 


A method has been proposed for the prediction 
of the CO, solubility in non-polar liquid hydro- 
carbons. In the gas phase the CO, solubility can 
be satisfactorily predicted using the Krichevsky— 
Ilinskaya equation combining with PR EOS + 
HV. For the liquid phase solubility of CO, in 
hydrocarbons, a similar method based provided 
accurate values of 05° are used. The method pro- 
posed does not require fitted interaction parame- 
ters or have problem of non-convergence. This 
study shows that accurate determination of the 
0x is critical for accurate prediction of solubility 
data at elevated pressures. The EOS generally 
over estimates the calculation of this parameter 
for most the systems studied. The 63 was opti- 
mised only to emphasise that the predictive accu- 
racy of the proposed model depends heavily on 
the accuracy in estimating this parameter. With 
0x estimated from an EOS with HV mixing rules, 





Deviation in 05 (%) 





Optimised 





0.0300 
0.0365 
0.055 
0.055 
0.010 
0.040 
0.009 
0.035 
0.001 
0.12 
0.12 
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the proposed model is predictive in nature, how- 
ever accurate method for the determination of i+ 
is necessary. With the use of optimised 05 values, 
the model becomes semi-predictive. 
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Abstract 


Fatty oil influence on the solubility of limonene in CO, was investigated under pressures 8—12 MPa at 313.2 K,a 
temperature typically applied in supercritical fluid extraction of essential oils. Solubility in CO, was measured using 
the dynamic method both for limonene and for the mixture of limonene and blackcurrant seed oil. In the whole range 
of pressures applied, the concentration of fatty oil in the vapour phase is negligible in comparison with the 
concentration of limonene. Limonene is distributed between the liquid phase rich in fatty oil and the vapour phase 
rich in CO,, and its equilibrium concentration in the latter decreases with the diminishing limonene-to-oil ratio in the 
saturator. There is a steep increase of the limonene partition coefficient with pressure between 8 and 10 MPa, near 
the critical pressure of the binary mixture of limonene and CO,. The observed behaviour of the three-component 
system was confirmed and explained by thermodynamic modelling. The thermodynamic model applied was the 
Soave—Redlich— Kwong cubic equation of state with either the one fluid linear van der Waals mixing rule or with the 
MHV?2 mixing rule. Appropriate conditions for an efficient supercritical fluid extraction of essential oils from seeds 
follow from the results obtained. Extraction pressure should be approx. 20% larger than the critical pressure of the 
essential oil + CO, binary mixture and rather tight packing of the ground seed in the bed should be applied. © 2001 


Elsevier Science B.V. All rights reserved. 
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supercritical fluid processing. For the design of an 
industrial plant, dependence of essential oil solu- 
bility on pressure and temperature is of interest. 
The frequently recommended extraction tempera- 
tures between 40 and 50°C and extraction pres- 
sures lower than 10 MPa are usually sufficient to 
sididdlaNasapees obtain the well-soluble volatile essential oils, while 

* Corresponding author. Tel.: + 420-2-20390234; fax: + the co-extraction of undesired higher molecular- 


420-2-20920661. ‘ yee 
E-mail address: sovova@icpf.cas.cz (H. Sovova). weight substances 1s limited. 


1. Introduction 


The present work relates to the extraction of 
essential oils with supercritical CO,, which is a 
promising field for the industrial application of 
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References to the data on the solubility of 
various essential oils components in CQO, are 
available in the survey by Reverchon [1]. How- 
ever, the solubility of essential oils extracted 
from air-dried aromatic plants is usually lower 
than their solubility measured in the pseudobi- 
nary system essential oil+CO,. This is con- 
nected with their low content in the material, 
ranging typically from tenths to several per cent. 
Identical solubilities of pure essential oil and of 
essential oil from plant material were reported 
only for the extraction of clove bud essential oil 
with an exceptional yield of 20.8% [2]. 

Phase behaviour of essential oils and other 
minor and trace components from plants, and in 
particular from leaves, has been investigated in 
connection with the analytical-scale supercritical 
fluid extraction (SCFE) for sample analysis. 
Hawthorne [3] pointed out that the phase equi- 
librium in such cases depends on distribution of 
the solute between the supercritical fluid and the 
sorptive sites on the vegetable matrix. Addition 
of an organic modifier to CO, greatly increases 
the extraction rates, indicating that matrix—so- 
lute interactions are more important for con- 
trolling the process than bulk solubility [4]. The 
‘matrix-effect’? was observed also by Goto et al. 
[5], who extracted essential oil from peppermint 
leaves. They simulated the dynamics of the su- 
percritical fluid extraction, using a model ac- 
counting for kinetics of the essential oil 
adsorption/desorption on the matrix. In a later 

















Fig. 1. Scheme of the apparatus for solubility measurement. 1, 
syringe pump; 2, pressure gauge; 3, saturator; 4, micrometer 
valve; 5, trap; 6, flow meter; 7, water bath; 8, cooling bath. 


work [6], the ‘matrix-effect’ was described by the 
BET adsorption isotherm. The isotherm enables 
a smooth transition from the extraction of con- 
centrated solute, where the equilibrium concen- 
tration in the solvent approaches the bulk 
solubility, to the extraction of trace amounts, 
where the equilibrium is determined solely by 
the solute interaction with the matrix. 

Extraction of essential oils from aromatic 
seeds is influenced by the presence of fatty oils 
(vegetable oils). The seeds contain tens of per 
cent of fatty oil as a store of energy for the 
future plant. Fatty oils are well miscible with 
essential oils, but their solubility in supercritical 
CO, is substantially lower than that of essential 
oils. A large decrease in the solubility of 
caraway essential oil extracted with CO, from 
seed was reported as a result of the ‘fatty oil 
effect’ [7]. 

In this work the effect of fatty oil was studied 
on a model system where the essential oil is 
represented by limonene, a common monoter- 
penic hydrocarbon. Its solubility in CO, was 
measured using the dynamic method both for 
the binary system and for the system containing 
fatty oil. Interpretation of experimental results 
would not be possible without a notion of the 
liquid-phase equilibrium composition of the 
ternary mixture, which however could not be 
obtained by the dynamic method. Phase equi- 
librium in the system was therefore analysed by 
means of thermodynamic modelling. 


2. Experimental 


2.1. Materials 


S(—)-limonene (97%) was purchased from 
Fluka Chemie AG, CO, (99.9%) was supplied 
by Technoplyn Linde. Fatty oil was extracted 
with CO, at 28 MPa and 40°C from blackcur- 
rant seeds supplied by a co-operative farm 
Chel¢ice. It was refined on a chromatographic 
column, so that it was composed only of triglyc- 
erides. Its effective carbon number 53.1 is close 
to the carbon number of triolein, 54. 





H. Sovova et al. /J. of Supercritical Fluids 20 (2001) 113-129 115 


Table | 
Solubility (mole fraction) of pure limonene in dense CO,: 
experimental data (Vjim,exp) and values calculated using Eq. (1) 


(Viim,caic) 





T (K) as (MPa) d (kg m~°) Viim.exp Viim,cale 





0.00144 
0.00123 
0.00190 
0.00216 
0.00292 
0.00320 
0.00396 


0.00135 
0.00150 
0.00183 
0.00218 
0.00266 
0.00313 
0.00367 
0.00418 
0.00169 
0.00201 
0.00313 
0.00366 
0.00404 
0.00524 
0.00674 


io ees 198.2 
Tod 210.7 
‘toa 232.0 
7.7 248.7 
7.88 266.0 
8.00 279.1 
8.10 291.2 
8.17 300.5 0.00431 
7.0 172.0 0.00174 
7.5 194.0 0.00183 
8.38 241.5 0.00285 
8.60 255.8 0.00345 
o.72 264.2 0.00383 
9.00 285.0 0.00485 
9.21 303.0 0.00594 
9.23 304.8 0.00658 0.00692 
9.30 311.1 0.00683 0.00761 
9.56 336.3 0.0108 0.0115 
9.59 339.4 0.0119 0.0121 





2.2. Experimental method 


Solubility measurement was performed in a 
flow-through apparatus (see Fig. 1). Carbon diox- 
ide was supplied to the equipment by a HPLC 
syringe pump. Saturator was a_ high-pressure 
column of i.d. 8 mm and volume of either 4 or 12 
cm’, filled with glass beads wetted with limonene 
or with a mixture of limonene and vegetable oil. 
The beads occupied more than half of the volume 
of the saturator. Glass wool and stainless steel 
frits placed at both ends of the saturator served as 
flow distributors and prevented liquid entrain- 
ment. The initial amount of limonene in the satu- 
rator ranged from 0.1 to 0.9 g, the maximum 
amount of oil was 0.53 g. Temperature was main- 
tained within + 0.1 K by a water bath. Pressure at 
the inlet of the saturator (at its bottom) was 
measured with accuracy + 0.01 MPa by a pressure 
gauge and controlled within + 0.05 MPa by the 
syringe pump. Pressurized parts of the equipment 
were connected by short capillaries of negligible 
volume. 


After 15 min of static extraction, CO, flow rate 
through the saturator was adjusted to 0.04—0.06 g 
min~'. Preliminary experiments were carried out 
at different flow rates using the saturator volume 
of 4 cm? to show that this flow is sufficiently slow 
to obtain full saturation of the solvent. The solu- 
tion from the saturator expanded in a heated 
micrometer valve to atmospheric pressure, and 
limonene separated from the gas was collected in 
a trap, a U-shaped glass tube equipped in the exit 
shoulder with a glass-wool filter in order to pre- 
vent entrainment of tiny limonene droplets. The 
trap was cooled in a bath containing a mixture of 
dry ice and ethanol. Gaseous CO, flowed from 
the trap to a wet flow meter before being vented. 

After a period of dynamic extraction the sol- 
vent flow was interrupted and the amount of 
limonene collected in the trap was determined by 
mass (typically 0.01—0.04 g, accuracy better than 
2%). The volume of passed CO, (typically 0.3—1.5 
litres, accuracy better than 1%) was converted to 
mass and limonene-to-CO, mass ratio was calcu- 
lated. The amount of limonene remaining in the 
saturator was calculated from mass balance. This 
procedure could be employed also when the mix- 
tures of limonene and fatty oil were measured, 
because the solubility of oil in CO, was negligible 
under the conditions applied, and the trap con- 
tained practically pure limonene. Finally, the eval- 
uated composition of the loaded solvent and the 
limonene-to-oil ratio in the saturator were ex- 
pressed as mole fractions, using in case of fatty oil 
the molecular weight of triolein. 

Experimental error of a single solubility mea- 
surement determined from repeated experiments 
was approx. 10%. The error is caused mainly by a 
small but variable hold-up of limonene condens- 
ing inside the micrometer valve. All results listed 
in Tables 1—5 are mean values from at least two 
or three measurements. 


2.3. Results 


The difference between the limonene concentra- 
tions in CO, measured for the binary and ternary 
systems is illustrated in Fig. 2, where the concen- 
tration is plotted against the average mass of 
limonene in the saturator during the sampling 
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Table 2 Table 5 
Limonene solubility (limonene-to-CO, mole ratio) at 8 MPa Dependence of limonene solubility (limonene-to-CO, mole 
and 313.2 K in presence of fatty oil ratio) on fatty oil at 12 MPa and 313.2 K 








Z\im/ Zoil Viim/ ¥co2 Zim! oil Yim Voo2 Z1im/ oil Viim/ Yco2 











0.23 0.0005 0.17 0.002 0.45 0.002 
0.45 0.0008 0.61 0.006 0.94 0.004 
0.83 0.0012 1.10 0.011 1.54 0.008 
1.41 0.0017 1.61 0.018 2.21 0.011 
2.21 0.0020 2.38 0.025 5.33 0.028 
2.96 0.0024 3.35 0.037 6.38 0.032 
3.86 0.0027 4.32 0.044 7.44 0.041 
4.95 0.0032 O=0.25 g O=0.15 g 

5.99 0.0034 S = 0.0106 S = 0.0053 

O = 0.09-0.31 g 
S = 0.002 

















pling was repeated in both experiments until prac- 
period. The figure depicts the results of two exper- tically all limonene was extracted. The horizontal 
iments performed at 8 MPa and 40°C, using the line indicates the experimental solubility of pure 
saturator volume of 4 cm*. Vapour phase sam- limonene in CO,, y = 0.0036. 


Table 3 
Dependence of limonene solubility (limonene-to-CO, mole ratio) on fatty oil at 9 MPa and 313.2 K 











7 7 eee i ee: ae nae ? | = A ly 
Ziim/Zoil Viim/ Vco2 “lim/ “oil J lim/ ¥co2 “lim/ <oil Viim/ ¥co2 








0.36 0.003 0.30 0.002 0.36 0.001 

0.63 0.005 0.64 0.003 0.86 0.002 

0.99 0.008 1.22 0.006 2.22 0.007 

1.35 0.012 1.69 0.009 3.54 0.011 

2.29 0.018 2.14 0.011 4.76 0.016 

2.47 0.018 2.75 0.015 6.19 0.021 

3.84 0.030 4.68 0.024 7.33 0.024 
6.17 0.033 9.64 0.035 
8.23 0.046 

O=0.27 g O=0.13 g O=0.07 g 

S = 0.0077 S = 0.0054 S = 0.0035 














Table 4 
Dependence of limonene solubility (imonene-to-CO, mole ratio) on fatty oil at 10 MPa and 313.2 K 











= |» ly) - ” ) ly) = = ’ , 
Z\im/Zoil Viim/ Yco2 Ztim/ Zoi Niim/ Vco2 Ziim/ Zoil Viim/ Voo2 














0.42 0.004 0.19 0.001 0.36 0.001 
0.92 0.010 0.34 0.002 1.12 0.007 
1.58 0.016 0.55 0.004 2.26 0.011 
225 0.025 1.15 0.009 3.72 0.021 
2.74 0.030 2.65 0.019 6.45 0.039 
kX 0.035 3.86 0.041 11.40 0.070 
3.86 0.041 

O=0.45 g O=0.13 g O=0.10 g 

S = 0.0106 S = 0.007 S = 0.006 
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For the binary system, the measured concentra- 
tion was initially above this line because the liquid 
phase filled the saturator up to the top and flowed 
to the trap together with the solvent. When the 
liquid level in the saturator dropped, the mea- 
sured concentration was equal to the solubility of 
limonene in CO). Finally, when the liquid phase 
in the saturator disappeared, the vapour phase 
was diluted by fresh CO, and the limonene con- 
centration dropped gradually to zero. In order to 
extend the exploitable number of measurements in 
One experimental run, the larger saturator volume 
of 12 cm? was used for these experiments. 

In the ternary system, the limonene solubility 
was initially equal to the solubility of pure 
limonene in CO,, and later, as the amount of 
limonene in the saturator diminished, it was de- 
creasing. Though the initial amount of limonene 
in the saturator was similar to that in the binary 
system, no liquid-phase overflow was observed, 
obviously because the density of the equilibrium 
fatty oil + limonene + CO, phase was higher than 
the density of the liquid phase in the binary 
system. 


2.3.1. Limonene + CO, 

The experiments were carried out at tempera- 
tures 313.2 and 323.2 K and at pressures mostly 
above the critical pressure of CO,, 7.4 MPa. The 
experimental solubilities listed in Table 1 were 
compared with literature data as follows. 





0,016 


0,012 + 


0,008 - 


y(limonene) 














0,1 0,15 0,2 
Limonene in saturator, g 


Fig. 2. Dependence of the limonene concentration (mole frac- 
tion) in the CO, exit stream on the amount of limonene in the 
saturator at 8 MPa and 313.2 K. Binary system: CL], exp. data, 
— , limonene solubility in CO,; ternary system: A exp. data 
for 0.21 g oil in the saturator. 
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The first step was to smooth and interpolate the 
solubility data using the equation derived by 
Adachi and Lu [8]. With its coefficients adjusted 
to our data the equation was found acceptable to 
fit the limonene solubility data of Iwai et al. [9], 
Di Giacomo et al. [10], Matos et al. [11] and 
Marteau et al. [12]. These data together cover a 
range of temperatures from 308.2 to 333.2 K, 
pressures from 3 to 10 MPa, and limonene mole 
fractions in CO, from 0.0005 to 0.023. Thus, as a 
second step we determined, minimizing the sum of 
squares of relative deviations between measured 
and calculated mole fractions, new coefficients of 
the equation from a combined set of the above 
literature data and our data from Table |. The 


result is 


C = 2-84 —0.00271d + O.00000 7050" orcas — 4409/T + 0.24) 
(1) 


where c is the limonene concentration in the su- 
percritical CO, (kg m~*), d is the pure CO, 
density (kg m~*) and T is the temperature (K). 
The CO, density was calculated according to the 
Altunin and Gadetskii equation of state [13], the 


mass fraction of limonene in carbon dioxide was 
calculated as w = c/d. The average absolute devia- 
tion of the calculated solubilities from our experi- 
mental results is less than 7%, 

Critical pressure of the limonene + CO, mixture 
at 313.2 K was estimated by Budich [14] to Pan = 
8.44 MPa. Limonene is fully miscible with CO, 
above this pressure, so that the concentrations 
measured in the binary system at 9 MPa, depicted 
in Fig. 3, indicate how the homogeneous mixture 
in the saturator is diluted with fresh CO>. The 
dilution data enabled us to describe by a simple 
model the solvent flow pattern in the saturator, 
or, more exactly, in the space between the closed 
valves at the saturator inlet and outlet (see Fig. 1). 

The model consists of an ideal mixer represent- 
ing the saturator, and a plug-flow section repre- 
senting the outlet part of the saturator together 
with the pressurized inlet part of the micrometer 
valve. Mass of the solvent in the mixer, M@,,, and 
in the plug-flow section, M,, are assumed to be 
constant during the dilution. Initial mass of 
limonene in the mixer and in the plug-flow section 
is Lio and Ly, respectively. Mass of the solvent 
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y(limonene) 














0,1 0,15 
Limonene in saturator, g 


Fig. 3. Dependence of the limonene concentration (mole frac- 
tion) in the CO, exit stream on the amount of limonene in the 
saturator at 9 MPa and 313.2 K. Binary system: A, exp. data, 
x average concentrations during the sampling intervals calcu- 
lated with the model of solvent flow, — , instantaneous 
concentration calculated with the model; ternary system: A, 
exp. data for 0.16 g oil in the saturator. 


passed through the saturator is Q, mass of the 
limonene passed from the mixer to the plug-flow 
section is E,,, and mass of the limonene passed 
from the micrometer valve to the trap is &,. 
According to the model, 


Em = Lmoll — exp(— Q/M,,)] (2) 

E,, = L,oQ/M, for O< M,, (3) 

E, = Loo + Lmol! — exp(— (Q — M,)/M)I for Q 
> M,. (4) 


Limonene-to-CO, ratio in the stream entering the 
trap, expressed in mass, is 


dE,/dQ = L,»/M, for O< My, 
dE,,/dQ = (Lino/M,,)exp( — (Q — M,)/M,,) for Q 
> M,,. (5) 


The initial mass of limonene in the whole section, 
Loo = Lmo + Lpo, is known. Mass of CO, in the 
section, M,= M,, + M,, was determined from the 
volume of CO, passed through the flow meter 
during depressurization of the section multiplied 
by the CO, density at ambient pressure. An equa- 
tion containing only two unknown quantities, L,o 
and M.,, is obtained by substituting these relations 
into Eq. (4): 


ED ss L5o ss (Leo aes Lo) 
[1 Fae exp( ae (Q ae M,,)/(M, SG M,))I for OQ 

> M,. (4a) 
The values Lo and M, were adjusted by fitting 
Eqs. (3) and (4a) to the measured amounts of 
limonene in the trap. Finally, volumes of the 
mixer, V,,, and of the plug-flow section, V, = 0.34 
cm*, were obtained from the known amounts M,,,, 
M,, and CO, density, d. This value of V,, was used 


in the evaluation of the experiments with ternary 
system. 


2.3.2. Limonene + oil + CO, 

The experiments were carried out at tempera- 
ture 313.2 K and at pressures 8, 9, 10 and 12 
MPa. Figs. 2 and 3 show that when oil is present 
the limonene solubility strongly depends on the 
content of limonene in the saturator and thus on 
the equilibrium liquid-phase composition which 
could not be measured by the dynamic method. 
Instead, we calculated from mass balance the 
equilibrium ratio of limonene to oil in the section 
represented by the mixer where phase equilibrium 
is supposed to be established. 

The ratio is a function of the amount of CO, 
passed through the saturator, Q. It follows from 
Eq. (5) that, except for the initial period O< 
M,=V,d, the composition of the micrometer 
valve exit stream at Q equals the vapour phase 
composition which was in the mixer at O— M,,. 
Further it follows from Eqs. (2) and (4) that the 
amount of limonene in the mixer at Q— M, was 
equal to L.,— E,(Q), the difference between the 
initial mass of limonene in the whole section and 
the mass of limonene in the trap at Q. Thus, from 
the amounts of limonene in the trap E,,, E,»...., 
E,, measured at Q,, Q3,..., Q,, only those sampled 
at Q;> M, were taken into account. The mean 
value of the equilibrium limonene-to-CO, ratio in 
the vapour phase for a sampling period was calcu- 
lated as (E,;,,—£,))/(Q;,;—Q,). The corre- 
sponding overall limonene-to-oil ratio in the 
saturator was calculated as [L.9 — (Ey; + £5: + 1)/2]/ 
O, where O is the mass of oil in the saturator. The 
mass ratios were converted to mole fractions, 
denoted in Tables 2—5 as yjin/V¥co, for the vapour 
phase and 2,,,,/Z,;, for the mixture in the saturator. 
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Dependence of Vjim/Vco, ON Ziim/Zco, WaS exam- 
ined for different amounts of oil in the saturator, 
O. No effect of the amount of oil was observed at 
8 MPa. Limonene solubility grew with increasing 
limonene-to-oil ratio and close to 2;,,/Z,; =5 it 
reached the limonene solubility in the binary system 
(see Table 2). The initial slope of the dependence, 
or the limonene partition coefficient at low 
limonene-to-oil ratio, was S = (\im/Vco,)/(Zim/ 
Zou) = 0.0020. ; 

Above P.,, at pressures 9-12 MPa, the limonene 
solubility dependence was more complex (see Ta- 
bles 3—5). For a constant amount of oil in the 
saturator, O, the limonene-to-CO, ratio in the 
vapour phase was directly proportional to the 
limonene-to-oil ratio, so that the data were charac- 
terized by a constant limonene partition coefficient 
S. With a growing amount of oil in the saturator 
the coefficient increased, and the maximum parti- 
tion coefficient was obtained just before the liquid 
phase started overflowing from the saturator. These 
experimental values are listed in the first columns 
of the tables. The maximum limonene partition 
coefficients are S= 0.0077 for 9 MPa, and S= 
0.0106 both for 10 and 12 MPa. 

The different behaviour at conditions below and 
above the critical point of the limonene + CO, 
mixture can be explained using the information 
obtained by modelling the phase behaviour of the 
three-component system. To adjust the parameters 
of the thermodynamic models employed, experi- 
mental data on the equilibrium composition of 
both vapour- and liquid phases in the constituent 
binary systems were necessary. From the data 
reported in the literature we selected those that 
closely agree with our experimental results on the 
vapour phase composition. For the system 
limonene + CO, it is the data of Iwai et al. [9]. For 
the system blackcurrant seed oil + CO, (solubility 
of this oil in CO, is reported elsewhere [15]) we used 
the data obtained by Klein and Schulz [16] for 
refined rape seed oil and by Bharath et al. [17] for 
sesame oil and for triolein. 


3. Thermodynamic modelling 


Cubic equations of state (CEOS) appear to be 
attractive to model the phase equilibria of essential 


oils and/or fatty oils + SC CO, systems. For exam- 
ple, Brandani et al. [18] correlated the solubility of 
pure limonene and citral in SC CO, using an EOS 
with density-dependent mixing rules. Then the 
equation was used to predict the phase diagram of 
the CO, + limonene system at temperatures from 
308.2 to 323.2 K and over the whole range of 
pressure and composition. Vieira de Melo et al. [19] 
applied Peng-Robinson (PR) CEOS with five dif- 
ferent mixing rules to calculate the vapour—liquid 
equilibria of CO, + limonene, and carried out a 
comparative analysis of their performance. PR 
CEOS combined with the one-fluid van der Waals 
mixing rule was also applied to calculate the phase 
equilibria of CO,+limonene in other works 
[9,20,21]. 

Mathias et al. [22] correlated by means of the PR 
EOS the solubility data of palm oil, a mixture of 
triglycerides, in high-pressure CO,. Geana and 
Steiner [23] also used the PR EOS with a conven- 
tional mixing rule involving two parameters to 
model the extraction of rapeseed oil with supercrit- 
ical CO,. Bottini et al. [24] applied the Group 
Contribution Equation of State (GC-EOS) to cal- 
culate the vapour—liquid, and liquid—liquid equi- 
libria of supercritical alkane—vegetable oil 
mixtures. Weber et al. [25] applied the Redlich— 
Kwong-—Aspen equation of state to calculate the 
vapour-—liquid equilibria of tristearin, tripalmitin, 
and triolein in CO, and propane. Coniglio et al. [26] 
correlated supercritical fluid—liquid equilibria for 
CO, and fish oil related compounds by means of 
the SRK EOS. Bharath et al. carried out an 
extensive experimental study and correlated the 
VLE of binary mixtures of fatty acids and their 
triglycerides in supercritical CO, using the PR 
CEOS [27]. 

In our study we use the Soave—Redlich— Kwong 
(SRK) CEOS with the one fluid van der Waals 
mixing rule to model the phase behaviour of the 
limonene + SC CO, and triolein + SC CO, binary 
systems: 


RT a 
v—b_ v(v+b) 
The one fluid van der Waals mixing rule is: 


Anix = yi xX; ay (7) 
| 





yr? (6) 
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Dynix = 2d (8) 


where a,,;, and b,,;, are, respectively, the mixture 
energy and co-volume parameter of the EOS. 
Usually, a geometric mean rule is applied to de- 
termine the cross-energy parameter a,j: 


= (a;;a aij ) (1 - kj) (7a) 


and either an arithmetic mean (conventional lin- 
ear rule) 
No 
Dinix = 2% bi: (8a) 
Pye 
or a quadratic rule for the cross-co-volume 
parameter 


b= (EPH) — 4) (8b) 


We apply the conventional linear mixing rule [Eq. 
(8a)], and thus only one binary interaction 
parameter has to be optimized to tune the EOS. 

To provide an additional basis for comparison 
and flexibility of the modelling, an alternative 
mixing rule for the mixture energy parameter is 
applied, namely the MHV2 model. In our case it 
is based on incorporating the modified UNIFAC 
[28] excess Gibbs energy model expression at zero 
pressure in the attractive term parameter a,,;, of 
the SRK CEOS [29,30]: 


N N 
12 se »y *) - toa “ y zai) 


i=] i=1 


aft ——e a nl “as (9) 


j=] ll 


In the above equation, z, is the phase composition 
in mole fractions, and the dimensionless parame- 
ters «,,;,and «are related to the mixture and pure 
compound parameters, 


a a:; 
= and a, =——, (9a) 


Amix = i a b RT 


mix 


~ -b, RT 


The recommended values of g, and g, are — 0.478 
and — 0.0047, respectively. 


3.1. Limonene + CO, 


The critical properties and acentric factor of 


d-limonene used in our calculations are those 
given in the Pure Component ae of the 
process simulator HYSYS [31], namely: 7. = 650 
K, P,=2.75 MPa, w = 0.3123. 

The phase equilibrium of the limonene + CO, 
system was modelled at three temperatures T= 
313.15 K, 323.15 K, and 333.15 K. Owing to a 
lack of space, a comparison between the measured 
and calculated mole fractions of CO, in the liquid 
and vapour phase as a function of pressure is 
presented here only for 7= 313.15 K (Fig. 4). 

The values of the temperature-dependent binary 
interaction parameter Kjim.co, [Eq. (7a)] were ob- 
tained via a nonlinear regression procedure ap- 
plied to the sets of experimental data measured by 
Iwai et al. [9], and are given in Table 6. The 
calculations demonstrated that the vapour—liquid 
equilibria of the limonene + CO, system is pre- 
dicted very well with a CEOS with a one fluid van 
der Waals conventional linear mixing rule for the 
mixture co-volume parameter. The results of the 
modelling represent adequately the experimental 
data [9] both in the liquid, and in the vapour 
phase. The maximum AARD of the SRK CEOS 
predictions from these data were calculated to be 
5.2% for the liquid phase, and 0.022% for the 


vapour phase. 
3.2. System CO, + triolein 


To apply the SRK CEOS to model the phase 
equilibrium of the CO, + triolein system, with 
either a one fluid van der Waals mixing rule or a 
mixing rule based on an excess Gibbs energy 
model, the critical parameters and acentric factor 
of triolein are required. For this high molecular 
weight compound this information is experimen- 
tally unattainable because of its negligible vapour 
pressure and thermolability. Thus, the parameters 
have to be theoretically estimated, and their val- 
ues taken as characteristic parameters (or consid- 
ered as hypothetical quantities) rather than 
properties with any physical meaning. The ap- 
proach adopted in this study is the following: 

Firstly the hypothetical normal boiling point 
and standard liquid density of triolein were esti- 
mated with the methods developed by Cholakov 
et al. [32]. Then, using this input information, 
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Fig. 4. Vapour-liquid equilibria calculations for the limonene + CO, system at T = 313.15 K as a function of pressure. Comparison 
with experimental data: — , MHV2 model, — : — - — , SRK CEOS; k,, = 0.1087; Ml, Iwai et al. [9]; A, this work. 


several sets of triolein critical properties were 
estimated with the in-built methods of the commer- 
cial process simulator HYSYS [31]. In order to 
choose the set(s) to be used in the subsequent 
modelling of the three component system, the 
influence of triolein parameters on the CO, + tri- 
olein phase equilibria calculations was assessed 
both for the case of the one fluid mixing rule (with 
Kco,,triolein = 9), and of the MHV2 model. For the 
former the best predictions were obtained with the 
following values for triolein critical parameters: 
T, = 810.03 K, P, = 1.09927 MPa. The value of the 
acentric factor was optimised by us and calculated 
to be w = 2.7171. In the paragraphs to follow this 
set will be referred to as set A. 

With the MHV2 mixing rule model the best 
predictions were obtained with a different set of 
critical parameters for triolein: 7, = 860.46 K, 
P, = 1.69942 MPa; and an optimized value for the 
acentric factor w = 0.7309. In all further discus- 
sions this set will be referred to as set B. Owing to 
a lack of space these preliminary phase equilibria 
calculations are not represented graphically here. 


Secondly, we studied the influence of the mixing 
rule employed in the thermodynamic model. We 
performed phase equilibrium calculations for the 


Table 6 
Binary interaction parameters for CO,+ limonene system 





Temperature k;; li, 
(K) 


Reference 





0.09 
0.09 
0.09 
0.11 


313.2 
323.2 
333.2 
313 


Iwai et al. [9] 


Akgiin et al. 
[20] 

0.11 

0.11 

0.094781 


323 
333 
313.2 


— 0.020437 


Vieira de Melo 
et al. [19] 
— 0.033065 
0.101861 —0.035184 
0.112756  —0.023589 
0.1087 - 
0.106976 _ 
0.106 - 


323.2 0.101823 
333.2 
343.2 
313.2 
sad. 


333.2 


This work 
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Fig. 5. Phase equilibria calculations for the system triolein + CO, at T= 313.15 K. Influence of the thermodynamic model mixing 


rule: — - — - — , SRK with conventional linear mixing rule. keo, triotein = 


— 0.0064, and triolein characteristic parameters set A; — 


, SRK with MHV2 mixing rule and triolein characteristic parameters set B; @, Bharath et al. [7], Ml, Klein and Schulz [9]. 


CO, + triolein system at T= 313.15 K and T= 


333.15 K. A comparison between the experimen- 
tally measured and calculated, with either the 
conventional mixing rule or with the MHV2 
model, mole fractions of CO, in the liquid and 
vapour phase, as a function of pressure, is pre- 
sented here only for T= 313.15 K (Fig. 5). 

The binary parameter ko triolein IN the conven- 
tional linear mixing rule was estimated via a 
nonlinear regression procedure using the experi- 
mental data of [16,17] on the equilibrium compo- 
sitions of the system vegetable fatty oil (or 
triolein) + CO, in the temperature range 313.2— 
333.2 K. At pressures lower than 10 MPa the 
agreement between the calculated and experimen- 
tally measured CO, mole fractions in the liquid 
phase is not good, while for the vapour phase it is 
excellent up to P= 15 MPa. The predictions for 
the liquid phase improve considerably in the pres- 
sure range 10-25 MPa, but deteriorate for the 
vapour phase. Moreover, at very high pressures a 
homogeneous phase is predicted only (not shown 
in Fig. 5). 

With the MHV2 mixing rule the predictions for 
the liquid phase composition are considerably bet- 


ter than for the vapour phase, for which the mole 
fraction of CO, is greatly underestimated. One 
possible reason behind this is that we describe the 
triolein applying the following functional groups: 
(CH2), (CH3), (CH=CH), and the ester group 
(COO), and we use parameters to describe the 
interactions between the corresponding groups 
and the CO, with values that have been deter- 
mined by Dahl et al. [30] from binary VLE data 
involving light esters only (methyl, ethyl, butyl 
acetate and ethyl butyrate). 

One way to improve the MHV2 model predic- 
tions is to define the ester chemical function as a 
new group, and to determine interactions between 
the new group and the rest of the groups [24,33]. 
In our study we adopted a different approach, 
which will be discussed in the following 
paragraph. 


3.3. Limonene + triolein + CO, 


The three component system limonene + tri- 
olein + CO, is a strongly nonlinear complex sys- 
tem. The modelling with good accuracy of the 
shape of its phase diagram throughout the entire 
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composition space presents many difficulties, 
some of which are directly connected with the 
robustness and convergence behaviour of the 
numerical methods and the two-phase flash 
algorithms applied. This question, however, is 
well beyond the scope of the present paper, and 
will not be highlighted here. It is only necessary to 
point out that in all cases studied the algorithms 
and numerical techniques which we employ 
demonstrated a _ steady and _ nonoscillatory 
convergence and practically did not exhibit a 
tendency towards a strong attraction to the trivial 
solution, even at conditions close to the critical 
point. Details of the organization of our novel 
strategy to predicting and calculating two- and 
three phase equilibria can be found elsewhere [34]. 

In this study we have modelled the phase 
behaviour of the three component system at 
T= 313.15 K and four different values of the 
pressure: P= 8, 9, 10 and 12 MPa. In the case 
when the SRK CEOS is applied with the 
conventional mixing rule, the interaction 
parameters for the constituent binaries 
CO, + limonene, and CO,+triolein are those 
determined as explained in the previous 
paragraphs. The binary interaction parameter for 
the triolein + limonene was accepted to be equal 
to zero, though this approximation is not 
completely adequate, since the two components 
are very diverse. The reason behind this is that 
experimental data for the triolein + limonene 
binary at any temperature is not available, and 
hence an optimised value of the Ajim trioten Cannot 
be estimated. In all cases studied we apply set A 
for the triolein characteristic parameters. 

To improve the predictions for the MHV2 
mixing rule the combinatorial part of the modified 
UNIFAC model was revised in order to account 
for such large size molecule as triolein. This idea 
was discussed originally by Smith et al. [35], who 
treated the exponent in the combinatorial term as 
a fitting parameter. 

The combinatorial expression used in modified 
UNIFAC is: 


Iny? = n( ee (10) 


x; 


we 


with modified volume fractions: 


(11) 


J 
In our work we also treated the exponent in the 
combinatorial term as an optimization parameter. 
The best predictions were obtained with an 
exponent value equal to 0.8675, which is very 
close to what Smith et al. [35] found for the 
Lyngby version of the UNIFAC. 

The vapour-liquid phase diagrams at P=8 
and 10 MPa, calculated with the conventional 
linear and with the MHV2 mixing rule, are shown 
in Figs. 6 and 7, respectively. The figures also 
show a comparison between the experimentally 
measured and the calculated by us dependence of 
the ratio of limonene versus CO, in the vapour 
phase to the limonene versus triolein ratio in the 
liquid phase. Graphs similar to those depicted in 
Fig. 7 were obtained at 9 and 12 MPa, but owing 
to a lack of space they are not given here. 


4. Discussion of the oil effect 


4.1. Solubility measurement 


The behaviour of the system with oil during the 
limonene solubility measurement is explained on 
schematic ternary diagrams. Fig. 8(a) shows the 
situation when limonene is fully miscible with 
CO,. Through the overall composition of the 
mixture in the saturator A passes the tie-line BC 
connecting the equilibrium vapour- and liquid- 
phase compositions B and C, respectively, which 
are positioned on the vapour-—liquid binodal. The 
available experimental equilibrium data are the 
limonene-to-oil ratio in the saturator, indicated by 
the line DE, and, disregarding the negligible con- 
tent of oil in the vapour phase, the vapour-phase 
composition B. As the amount of CO, in the 
saturator is unknown, the point A cannot be 
located. The liquid-phase composition C is also 
not known. The variables Zjim/Zoi 294 Viim/Yco, in 
Tables 3—5 correspond to the line DE and to the 
point B, respectively. The more liquid phase con- 
tains the saturator, the closer is the point F to the 
equilibrium point C, and the closer is the overall 
mole ratio z,,,/z,, in the saturator to the equi- 
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librium liquid-phase ratio Xjim/Xou- The change of 
the overall composition from A to A’ results into 
a shift of the equilibrium liquid-phase composition 
to C’ and into an increase of the equilibrium 
vapour-phase limonene concentration from B to B’. 
In the experiments with the largest amounts of oil 
in the saturator, close to an overflow of the liquid 
phase, the points F and C were almost identical. 
These experimental points could therefore be com- 
pared with the results of the thermodynamic mod- 
elling. 

The situation at pressure P = 8 MPa is shown in 
Fig. 8(b). Limonene is only partially miscible with 
CO, and the single vapour phase zone is limited to 
a very small area at the top of the ternary diagram 
D. The tie-lines BC starting from points B are 
therefore overlapping the constant limonene-to-oil 
lines DE starting from point D. Point F, for which 
the limonene-to-oil ratio is known, is practically 
identical with the liquid-phase equilibrium point C 
no matter what the amount of the liquid in the 
saturator is, and z,,,,/z,;, is practically equal to 
Xtim/*Xoi- Independence of the vapour-phase compo- 
sition on the amount of liquid phase in the satura- 
tor was proved experimentally. Thus the data in 
Table 2 could be calculated as average values from 
experimental points measured with different 
amounts of oil in the saturator. 











limonene 


a) 


4.2. Consequences for the extraction from seed 


The solvent flow pattern during the solubility 
measurement was close to that of an ideal mixer, 
because the saturator was short and the flow was 
slow and interrupted. The conditions during super- 
critical extraction of essential oils are different. The 
fixed bed of the ground plant material in the 
extractor is much longer, the solvent flows faster 
and rather in a plug flow. Early after entering the 
bed the solvent becomes saturated with the solute, 
and then it flows through the bed without any mass 
transfer. When the material close to the solvent 
inlet becomes exhausted, the region of mass trans- 
fer shifts, until it reaches the end of the bed. Up to 
this moment the exit stream concentration remains 
constant and in equilibrium with the initial compo- 
sition of the extracted material; extraction rate is 
constant. 

The initial exit stream concentration is usually 
decisive for the process economics, as the largest 
part of solute is extracted in the first extraction 
period. The concentration of essential oil in the 
solvent should be maximized, keeping at the same 
time as low as possible the vapour-phase concentra- 
tion of other plant substances, such as fatty oil in 
the case of extraction from seed. Let us assume for 
the sake of simplicity that the properties of 





limonene 


Fig. 8. The equilibrium concentrations in the saturator at 313.2 K: (a) 9-12 MPa, (b) 8 MPa. 
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essential oil are identical with those of limonene. 
Content of essential oil in seed is typically 1—3%; 
content of fatty oil is 10-40%. The ratio of es- 
sential oil to fatty oil ranges from 0.025 to 0.3 by 
mass and from 0.16 to 2, when expressed in 
moles. The maximum essential oil concentration 
in carbon dioxide will be achieved at pressure 10 
MPa = 1.2P.,. Further increase in pressure would 
not affect the essential oil solubility; it could 
however enable co-extraction of the less soluble 
plant substances. 

Seed composition at the beginning of dynamic 
extraction differs from the original seed composi- 
tion, because a part of the essential oil dissolves 
in CO, while the extractor is pressurized. We 
assume that mass transfer takes place uniformly 
in the whole bed and that phase equilibrium is 
established before the solvent starts flowing from 
the extractor. (This is possible if a greater part of 
cell walls has been broken during the seed grind- 
ing.) If essential oil is substituted for limonene, 
Fig. 8(a) illustrates the situation in the extractor. 
The line DE represents the ratio of essential oil 
to fatty oil in the original seed. Point A repre- 
sents the overall mixture composition, point C 
— the composition of the liquid phase in the 
seed and point B — the equilibrium vapour- 
phase composition, decisive for the extraction 
rate. The more essential oil the overall mixture in 
the extractor contains, the higher is the rate. As 
the figure shows, for a seed of a given composi- 
tion such increase can be achieved by a shift 
from A to A’, which means by packing the bed 
of ground seed more tightly and leaving less 
space for CO, . 

The recommended extraction conditions are 
therefore a rather tight packing of the ground 
seed, and extraction pressure 1.2P.,, where P,, is 
the critical pressure of essential oil + CO, mix- 
ture at the extraction temperature. These condi- 
tions will be examined experimentally in a 
further work. 


5. Conclusion 


Our study has shown that fatty oil decreases 
the solubility of essential oils extracted with su- 


percritical carbon dioxide from seeds, compared 
to their solubility in the binary mixture with 
CO,. 

The solubility of pure limonene in supercritical 
CO, measured at 313.2 and 323.2 K agrees well 
with selected literature data. Decrease of the 
limonene solubility after addition of blackcurrant 
seed oil to the system was determined experimen- 
tally at 313.2 K and 8—12 MPa. 

Phase equilibria for the ternary system and for 
its binary constituents were simulated using the 
Soave—Redlich—Kwong equation of state with 
either the van der Waals classical mixing rules or 
the MHV2 model for the mixture energy 
parameter. The comparison between the predic- 
tions obtained with the two mixing rules illus- 
trates that the limonene—CO, interaction 
parameter in Eq. (7a) plays a more significant 
role than the choice of the one fluid van der 
Waals classical mixing rule over the MHV2 
model. 

The experiments combined with the thermody- 
namic modelling of the ternary system have 
demonstrated the difference in the phase _be- 
haviour below and above the critical pressure of 
the binary mixture limonene + CO,, P.,, and its 
impact on the supercritical fluid extraction. The 
limonene solubility in CO, decreases with dimin- 
ishing limonene-to-fatty oil ratio in the mixture. 
Above the critical pressure P., it decreases also 
with diminishing liquid-to-vapour ratio, while be- 
low P., it is independent of this quantity. The 
effect of pressure on the solubility is most pro- 
nounced in the vicinity of P.,, but above 1.2P., it 
becomes negligible. The recommended conditions 
for the extraction of essential oils from seed with 
CO,, based on these results, are extraction pres- 
sure near to 1.2P., and a rather tight packing of 
the bed of ground seed. 
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Abstract 


Solubilities of the fat-soluble vitamins K,, D; and D, in CO, and propane were measured at temperatures 30, 40, 
60 and 80°C and over a pressure range of 80—300 bar for CO, and 10—110 bar for propane using a static-analytic 
method. In order to verify weather, the solubility isotherm crosses the three-phase S—L—V line, solid—liquid (S—L) 
phase transitions of vitamins K;, D; and D, under the pressure of gas were determined with modified capillary 
method. For all vitamins under the pressure of CO,, propane and dimethylether, a negative dP/dT slope of the 
S—L-—V curve was found. For K; under the pressure of inert gas (argon and nitrogen), the S—L—V-curve had a 
positive dP/dT slope. The solubilities of vitamins D in dense CO, under the conditions investigated were in the range 
0.04 x 10~3—1.45 x 10~* mole fraction. For vitamin K, the solubility in CO, was higher and was in the range 
0.16 x 10~+—4.07 x 10~ ? mole fraction. The solubilities of vitamins D in propane were approximately up to 10 times 
higher as in CO, and in the range from 0.15 x 10~* to 12.40 x 10~* mole fraction. Opposite, the solubility of vitamin 
K,; in propane was lower than in CO, and was in the range from 0.03 x 10~+ to 2.99 x 10~? mole fraction. The 
experimental S-L—V data were correlated using Peng-Robinson equation of state and binary iteration parameters 
obtained from best fit were used for solubility calculation. The fit of S-L—V lines showed good results in the case of 
vitamin K, and the order of magnitude of the solid solubilities were well represented. © 2001 Elsevier Science B.V. 


All rights reserved. 


Keywords: Phase equilibria; Solubility; S-L—V line; Fat soluble vitamins; Dense gases 





Mostly, reported conventional methods for con- 
centration of vitamins from plant materials are 
difficult, costly or inefficient. Their industrial 
commercialisation is usually limited due to low 
yields concerning to costs. Supercritical fluid ex- 
traction (SFE) presents an alternative to conven- 
Bik cele OE tional separation processes. For the design of a 

* Corresponding author. Tel.: + 386-222-9446]; fax: + 386- separation process, the knowledge on phase equi- 


225-16750. i ; : At Bega 
E-mail address: zeljko.knez@uni-mb.si (Z. Knez). libria, which includes the equilibrium solubilities 


1. Introduction 


At the moment, fat-soluble vitamins are mostly 
produced synthetically. However, their isolation 
from natural sources is of increasing interest [1]. 


0896-8446/01/$ - see front matter © 2001 Elsevier Science B.V. All rights reserved. 
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Fig. 1. Equipment for determination of melting point under 
the pressure of gas with the Modified capillary method: (1) 
glass capillary; (2) high pressure cell; (3) lamp; (4) video 
camera; (5) computer; (6) safety valve; (7) high pressure pump. 
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Fig. 2. S-L—V phase equilibria for binary system D,—gas 
(CO,, propane or dimethylether). 
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Fig. 3. S-L-—V phase equilibria for binary system D,-gas 
(CO,, propane or dimethylether). 


of a solute in dense gas and solid—liquid phase 
transition of a substance under the pressure of 
gas, is essential. 

In general, four types of S-L—V phase be- 
haviour are known. When the SCF is relatively 
high soluble in the molten substance, the S-L—V 
curve has a negative dP/dT slope. In the second 
type, where the S—L—V curve has a positive dP/ 
dT slope, the SCF is only slightly soluble in the 
molten heavy component, and therefore, the in- 
crease in hydrostatic pressure raises the melting 
temperature. The third type, where both compet- 
ing effects of increasing solubility of the SCF in 
heavy component and increase in hydrostatic 
pressure are influencing the behaviour, the three- 
phase S—L-—V curve shows a temperature mini- 
mum. In the fourth type, which was first observed 
recently by Weidner et al. [2] for the system 
polyethyleneglycol-CO, and Weber et al. [3] for 
the systems tripalmitin—CO, and tristearin—CO,, 
the P—T projection of the S-L—V curve shows a 
temperature maximum and a_ temperature 
minimum. 

For the tocopherols (vitamin E) [1,4—9], vita- 
mins K, D, A [1] and B-carotene (provitamin A) 
[1,10—16] few solubility data in near- and super- 
critical CO, could be found in the literature. 
However, there were no data on the solubilities of 
vitamins K and D in propane and the S—L phase 
transitions of vitamins K;, D; and D, under the 
pressure of gas. 

In this work, the equilibrium solubilities of 
vitamins K,, D,; and D, were measured in CO, 
and propane at temperatures 30, 40, 60 and 80°C. 
In the case of CO,, isotherms were determined 
over a pressure range from 80 to 300 bar, whereas 
for propane experiments were carried out in the 
pressure range from 10 to 110 bar. 

Vitamins K,, D; and D, are solids at atmo- 
spheric conditions with relatively low melting 
points in the range from 86 to 118.5°C. Due to 
the fact, that melting point depression of a sub- 
stance can occur under the pressure of gas, it was 
verified weather the solid—liquid transition of vi- 
tamins occur along solubility isotherms and how 
the course of solubility isotherm is changed in this 


case. 
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Fig. 4. S-L—V phase equilibria for binary system K,—gas 
(CO,, propane, dimethylether, argon or nitrogen). 


The experimental S-L—V data were correlated 
using Peng-Robinson equation of state and binary 
iteration parameters obtained from best fit were 
used for solubility calculation. 


2. Experimental 


2.1. Determination of melting point under 
pressure 


The three phase S—L—V lines for vitamins (K3, 
D; and D,) under the pressure of gas (CO, 
propane, dimethylether, argon and nitrogen) were 
determined using modified capillary method (Fig. 
1). A detailed description of the experimental 
system and procedure can be found in literature 
[2,17]. The apparatus consists of the optical cell, 
which is designed for a pressure of 300 bar and 
temperature of 100°C. The volume of the cell is 14 
ml. The observation windows are made from sap- 
phire and are fixed to the cell with screws. The cell 
is equipped with three additional openings for 
introducing and emptying the gas and introducing 
a thermocouple. The thermocouple was calibrated 
by pure substances with known melting points 
(o-nitrophenol T,,=45.0°C, azobenzol T,, = 
68.0°C, bis-diphenyl ketone 7,, = 95.0°C, acetani- 
lide 7,,=115.0°C). The pressurised gas was 
introduced by a high pressure pump. The pressure 
was measured by electronic pressure gauge (Di- 


gibar PE 500 Hottinger-Baldwin to +0.1%) and 
the cell was electrically thermostated by a heating 
jacket to +0.5°C. The substance was filled in the 
glass capillary and put together with the thermo- 
couple into the high pressure cell. The gas was 
introduced to the cell and later heated. The melting 
of the substance in the capillary could be observed. 
Temperature and pressure were registered at the 
beginning and at the end of melting. In screening 
experiments it was checked whether the observed 
melting temperature range depended on the mode 
of operation (pressure/temperature increase: ‘up- 
ward strategy’ or pressure/temperature decrease: 
‘downward strategy’). No hysteresis was found 
within experimental accuracy ( + 0.2°C/ + 0.5 bar). 


2.2. Determination of the equilibrium solubilities 


For the solubility measurements a static analytic 
method was used. A detailed description could be 
found in literature [10,18]. The 120 ml equilibrium 
cell was loaded with a sufficient amount of the 
substance. The CO, from the supply tank was 
cooled to a liquid state and compressed into the 
equilibrium cell by a high pressure pump. The 
content of the autoclave was mixed with a magnetic 
stirrer under constant operating conditions (tem- 
perature and pressure) until the equilibrium was 
reached. After 1 h of phase separation, a sample 
(~ 0.2 ml) was taken through the sampling valve 
into a trap with a solvent (ethanol) where the 
components were solubilised. The trap was im- 
mersed in a sub-zero ethylene glycol-water bath. 
The expanding valve and lines were purged with 
~ 5—7 ml of solvent (ethanol). The amount of CO, 
released was measured with a disposal of water in 
graduated cylinder. The change in the pressure 
observed when taking the samples was from 0.5 to 
1 bar. A temperature change was not detected. 
Since the quantity of the sample was sufficiently 
small compared to the volume of the equilibrium 
cell, further experiments could be done. 

The concentration of the solute in ethanol was 
determined by UV spectrophotometry. The ab- 
sorbance was measured (accurate to + 0.001) at the 
absorption maximum of 265 nm for vitamins D, 
and D, and at 251 nm for K,. 
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2.3. Materials 


Vitamins were obtained from Aldrich (Deisen- 
hofen, D): vitamin D, of 99% purity (Cat. No. 
V380-8), D, of 98% purity (Cat. No. 27, 158-6) 
and K, of 98% purity (Cat. No. M5,740-5) and 
were used without further purification. CO, 
(99.995%), propane (99.5%) and dimethylether 
(99.9%) were obtained from Linde Plin (Celje, SI). 
Argon (99.998%) and nitrogen (99.996%) were 
obtained from MESSER (MG-Ruse, SI). 


Table 1 


3. Results and discussion 


3.1. Melting points under pressure 


The S—L Phase transitions of vitamins D,, D, 
and K, under pressure of gas (CO,, propane and 
dimethylether) are presented in Figs. 2—4 and 
Tables 1—3. For the comparison, S—L phase tran- 
sition of vitamin K, under pressure of inert gas 
(nitrogen and argon) was measured (Fig. 4, Table 
4). In the figures and tables, the mean value 


S-L phase transitions of vitamins K;,, D; and D, under pressure of CO, 





K; 
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Table 2 


S-L phase transitions of vitamins K,;, D, and D, under pressure of propane 





K; 




















Table 3 


S-L phase transitions of vitamins K;, D, and D, under pressure of dimethylether 








K, 








P (bar) 








8.0 
16.3 








between the upper and lower melting tempera- 
tures is presented. 

The melting point of vitamin D, (Fig. 2) under 
pressure of CO, (Table 1) decreases from 118.5°C 
at atmospheric conditions to 66.2°C (+4.6°) in 
the pressure range of 180—270 bar. The melting 
point depression of vitamin D, is even more pro- 
nounced under the pressure of propane (Fig. 2, 
Table 2) where it decreases to 52°C at 18 bar and 
remains constant with further raise of pressure. 

The melting point of vitamin D, (Fig. 3, Table 
1) under pressure of CO, decreases with the in- 
crease of pressure from 86°C at atmospheric con- 
ditions to 39°C at 130 bar and with a further 
increase of pressure only a very small increase in 
melting point could be observed (to 42°C at 300 
bar). Similarly to vitamin D,, vitamin D; melting 
point depression is even more pronounced under 
the pressure of propane (Fig. 3, Table 2) where it 
decreases to 35°C at 8.5 bar. Melting point de- 
pression of vitamin D, is the highest under the 
pressure of dimethylether (Fig. 3, Table 3), where 
it decreases to 16.6°C at pressure 3 bar. 


The melting point of vitamin K, (Fig. 4) under 
pressure of CO, (Table 1) decreases with pressure 
from 106°C at atmospheric conditions to 79°C at 
155 bar and with a further increase of pressure 
only a small decrease in melting point could be 
observed (74°C at 317 bar). Opposite to the case 


Table 4 
S—L phase transition of vitamin K, under pressure of nitrogen 


and argon 





Nitrogen Argon 





T (°C) T (°C) 





104.3 ’ 102.7 
104.6 103.2 
104.6 103.8 
105.3 104.0 
105.3 104.9 
106.2 105.1 
107.0 105.5 
107.9 105.5 
108.9 105.7 
109.9 105.6 
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ig. 5. Solubility of vitamin D, in: (a) CO,; and (b) propane. 
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Fig. 6. Solubility of vitamin D, in CO,: comparison with the 
literature data from Johannsen and Brunner [1]. 


of vitamins D, the melting point depression of 
vitamin K, is less pronounced in the presence of 
propane (Table 2). The melting point decreases 











A 
a a4 


O 


Fe] Oe me 
Bs ot. Me aS SS 


0 50 100 
(b) P (bar) 








Fig. 7. Solubility of vitamin D, in: (a) CO,; and (b) propane. 


1000 








g D; /kg propane 








p (kg/m3) 


Fig. 8. Solubility of vitamin D, in propane: g D3 vs. solvent 
density plot. 


with the increase of pressure to the minimum of 
94.9°C at 39 bar and with further increase of 
pressure it slowly increases to 98°C at 265 bar. 
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Fig. 9. Solubility of vitamin D, in CO,: comparison with the 
literature data from Johannsen and Brunner [1]. 
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Fig. 10. Solubility of vitamin K, in: (a) CO,; and (b) propane. 


Melting point depression of vitamin K, is the 
highest under the pressure of dimethylether (Fig. 


4, Table 3), where it decreases to 68.1°C at pres- 
sure 16.3 bar. S—L phase transition of vitamin K, 
under pressure of inert gas (nitrogen and argon) is 
presented on Fig. 4 and Table 4 and a raise of 
melting point with the increase of pressure can be 
observed due to low solubility of gas in vitamin 
K;. 

It can be concluded, that for the vitamins D,, 
D, and K, under pressure of gas (CO,, propane 
and dimethylether) in all cases the melting point 
temperature depression when increasing the pres- 
sure can be observed and is relatively pronounced 
at lower pressures (the S-L—V curve has a nega- 
tive dP/dT slope). At higher pressure the melting 
point temperature does not change any more and 
remains practically constant or rises again at 
higher pressures. 


3.2. Equilibrium solubilities 


The solubilities of vitamins D,, D; and K, in 
dense CO, and propane are presented on Figs. 
5—11 and Tables 5—7. Each data point represents 
the average of a least two measurements and the 
relative standard deviations ranged from 0.2 to 
15.6%. 

The solubility of vitamin D, in CO, (Fig. 5a, 
Table 5) is in the range from 0.043 x 10~? at 
80°C and 110 bar to 0.516 x 10~* mole fraction 
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Fig. 11. Solubility of vitamin K, in CO,: comparison with the 
literature data from Johannsen and Brunner [I]. 
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Table 5 
Solubility of vitamin D, in CO, and propane 








Propane 








p CO, (kg/m?) yx 10° g D,/kg CO, P (bar) p prop. (kg/m?) g D,/kg prop. 











701.4 0.053 0.475 8.9 4.847 
715.7 0.058 0.526 ' 19.3 18.984 
747.3 0.071 0.639 487.3 18.717 
812.2 0.126 1.140 492.1 22.869 
838.2 0.141 1.267 498.2 21.630 
877.9 0.242 2.185 502.1 21.234 
914.8 0.225 2.027 505.6 21.799 
931.4 0.203 1.831 508.5 20.893 
942.8 0.221 1.996 

945.5 0.246 2.215 


318.1 0.083 0.747 
647.0 0.094 0.848 
731.3 0.114 1.031 
793.9 0.134 LZ 
827.1 0.129 1.163 
857.7 0.167 1.503 
881.1 0.240 2.162 
883.8 0.220 1.984 
895.6 0.248 2.240 
910.5 0.466 4.205 
912.2 0.516 4.658 


207.9 0.173 1.557 
302.3 0.748 6.750 
493.7 0.825 7.445 
500.3 0.944 8.516 
506.7 1.325 11.960 
513.0 1.176 10.611 
525.0 0.912 8.233 
623.1 0.949 8.560 
693.7 0.723 6.518 
Totes 0.581 5.242 
792.2 0.529 4.771 
822.1 0.233 2.101 


259.2 0.043 0.384 
372.3 0.189 1.703 
469.1 0.276 2.484 
560.5 0.232 2.094 
621.7 0.305 2.750 
687.0 0.156 1.410 
726.5 0.109 0.984 
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at 40°C and 302 bar. Isotherms at 30 and 40°C and are in the range from 0.173 x 10~? at 84 bar 
present the solubility of a solid vitamin D, in CO, over maximum 1.33 x 10~* mole fraction at 130 
while 80°C isotherm presents the solubility of a bar to 0.23 x 10~* mole fraction at 289 bar. 
liquid vitamin in CO,. A course of isotherm at Unusual course of 60°C isotherm is probably a 
60°C is unusual and solubilities are approximately consequence of the fact that the isotherm lies very 
up to 10 times higher as at other temperatures close to S—L—V line (Fig. 2) and a S—L phase 


Table 6 
Solubility of vitamin D, in CO, and propane 





CO, Propane 





P (bar) p CO, (kg/m*) yx10? sg D,/kg CO, P (bar) p prop. (kg/m?) yx 10° g D;/kg prop. 





30°C 
9.0 18.9 1.508 13.201 
37.5 491.4 0.684 5.984 
54.3 496.1 1.625 14.243 
85.0 504.0 0.573 5.009 
102.5 508.0 1.662 14.553 
126.0 512.7 2.596 22.758 
144.0 515.9 1943 69.729 
40°C 
13.0 28.0 1.412 12.364 
47.0 478.6 1.884 16.499 
62.5 483.7 0.680 5.950 
76.0 487.9 0.758 6.635 
91.5 492.5 1.632 14.293 
94.0 493.2 1.612 14.120 
106.5 496.5 3.957 34.737 
132.8 502.9 3.422 30.017 
141.0 504.7 5.418 47.626 


60°C 
16.8 34.8 0.650 5.685 
20.5 426.5 6.855 60.341 
24.0 428.7 8.638 76.173 
33.0 434.2 6.909 60.821 
38.0 437.1 7.447 65.592 
46.0 441.6 10.898 96.320 
50.5 444.0 8.185 72.148 
66.0 451.9 10.660 94.191 
106.5 469.0 10.887 96.222 
115.5 472.1 12.398 109.745 


80°C 
24.8 $3.5 0.145 1.266 
34.0 347.7 0.178 1.559 
50.5 395.4 5.009 44.048 
69.5 412.8 7.206 63.454 
74.5 416.9 7.062 62.177 
78.5 420.0 6.056 53.260 
86.5 425.8 8.306 73.220 
102.0 435.4 7.343 64.667 
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Table 7 
Solubility of vitamin K, in CO, and propane 





CO, Propane 





P (bar) p CO, (kg/m?) yx 10° g K,;/kg CO, P (bar) p prop. (kg/m?) yx 10° g K,/kg prop. 





303.1 
575.0 
784.0 
838.3 
882.5 
914.6 


197.5 
260.9 
345.0 
369.2 
600.8 
718.4 
738.5 
782.1 
832.8 


164.5 
311.2 
450.8 
600.2 
682.9 
745.1 


0.466 
1.152 
2.067 
3.423 
4.068 
3.809 


2.44 
3.43 
kB 
3.04 
1.62 
3.12 
3.38 
3.15 
sat 


0.158 
0.484 
0.740 
1.325 
2.040 
2.688 


40°C 


1.824 
4.514 
8.105 
13.439 
15.984 
14.962 


11.5 
18.5 
42.0 
60.0 
83.0 
103.0 


60°C 


9.575 
13.479 
12.631 
11.939 

6.351 
12.241 
13.263 
12.374 
12.841 


0.617 
1.896 
2.898 
5.969 
8.001 
10.547 


18.0 
ft 
44.0 
65.5 
83.0 


23.9 
468.2 
476.8 
482.9 
490.0 
495.6 


38.6 
429.6 
440.5 
451.7 
459.7 
469.0 


0.025 
0.113 
0.202 
0.229 
0.191 
0.159 


0.045 
0.377 
1.094 
1.263 
1.155 
0.788 


0.097 
0.441 
0.789 
0.898 
0.749 
0.622 


0.177 
1.476 
4.287 
4.950 
4.525 
3.087 





Table 8 


Physico-chemical properties of vitamins K,, D, and D, 





Property 


K; 


D, 


D; 





Molecular weight 
Boiling point 
Critical temperature 
Critical pressure 
Critical volume 
Acentric factor 
Density 


Melting point 
Fusion enthalpy 
Fusion entrophy 


M,, (g/mol) 
T, (°C) 

T. (°C) 

P. (bar) 

V. (m3/kmol) 
a) 

p (g/cm?) 


(°C) 
AH (J/mol) 
AS*™s (J/mol K) 


FORMULA 
JOBACK* 
LYDERSEN®* 
LYDERSEN?* 
LYDERSEN?* 
DEFINITION? 


IMMIRZI AND PERINI? 


Measured‘ 
Measured‘ 
Measured‘ 
Measured4 


172.18 


365.05 
624.81 


36.079 
0.4880 
0.6385 
1.29213 
1.285942 

106.3 
21 488.06 
56.63 


396.66 
715.76 
938.27 
11.870 
1.3440 
1.0569 
1.0396 
0.96996 
117.9 
37 322.14 
95.44 


384.64 
689.16 
905.18 
12.136 
1.3130 
1.0808 
1.035 
0.97248 
81.1 
17 490.35 
49.37 





* ASPEN software [23]. 


> Calculated with method of IMMIRZI IN PERINI [24]. 


© Measured with helium pycnometer. 
4 Measured with DSC (Thermoanalytical system Mettler TA 3000). 
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Table 9 


Binary interaction parameters k,; and c,, for the system vitamin K,—dense gas (CO,, propane) and average absolute relative deviation 


AARD of calculated from experimental values 








System k;; Nij 


AARD (%)? 











S-L-V 








0.0121 
0.1109 


K,-CO, 
K ,—-propane 


1.02 
0.22 














® AARD(%) = 100. pv, Mate. =Yerol 
N Yexp 
transition of vitamin influences the solubility mea- 
surements. Similar unusual course of the solubility 
isotherm was already observed by Knez and Steiner 
[18] for the system capsaicin—CQ,. 

The solubility of vitamin D, in propane (Fig. 
5(b), Table 5) is higher than in CO,. It is almost 
independent from the temperature and increases 
with pressure from 0.22 x 10~? mole fraction at 9.8 
bar and 40°C to 2.79 x 10~3 mole fraction at 35 
bar and 60°C and remains in the range of 2.5 x 
10-3 mole fraction for all temperatures with a 
further increase of pressure to 110 bar. Isotherms 
at 30 and 40°C present the solubility of a solid 
vitamin while 60 and 80°C isotherms present the 
solubility of a liquid vitamin in propane (Fig. 2). 
No differences in course of solubility isotherms can 
be observed. Fig. 6 presents the comparison of 
experimental solubilities of vitamin D, in CO, with 
the literature data [1]. A strong disagreement be- 
tween two sets of data can be observed especially 
at 60°C, what could be a consequence of S—L phase 
transition of vitamin. 

The solubility of vitamin D; in CO, (Fig. 7(a), 
Table 6) is in the range from 0.24 x 10~? mole 
fraction at 80°C and 98 bar and 1.45 x 10~? mole 
fraction at 80°C and 294 bar. 60 and 80°C 
isotherms present the solubility of liquid vitamin 
while 40°C isotherm lies on the S-L—V phase line. 
The course of 40°C solubility isotherm is not 
influenced by the S—L transition of vitamin D,. 

The solubility of vitamin D, in propane (Fig. 
7(b), Table 6) is in the range from 0.15 x 10~? mole 
fraction at 80°C and 25 bar to 12.4 x 10~? mole 
fraction at 60°C and 116 bar. All solubility 
isotherms present the solubility of a liquid vitamin 
in propane. Solubility at 30 and 40°C is indepen- 


dent from temperature and is constant up to 100 
bar (0.6 x 10~3—1.9 x 10>? mole fraction at 40°C 
and 47 bar) and increases with further increase of 
pressure. When solubilities are presented on y 
versus P plot, the solubilities at 80°C are lower as 
at 60°C. However, when solubilities are presented 
as weight fractions vs. density of solvent, the 
solubility increases with increasing solvent density 
at constant temperature and with increasing tem- 
perature at constant density within the propane 
density range from 400 to 520 kg/m? (Fig. 8). 

Fig. 9 presents the comparison of experimental 
solubilities of vitamin D, in CO, with the literature 
data [1], and good agreement can be observed. The 
difference calculated as |p! — ythis work|/y this work ig 
in the range 5.9—30.4%. 

The solubility of vitamin K, in CO, (Fig. 10(a), 
Table 7) is in the range from 0.16 x 10~* mole 
fraction at 80°C and 82 bar to 4.1 x 10~* mole 


500 7 
450 - 
400 + 
350 - 
—<>300 4 
£250 4 
00 4 
150 4 
100 - 
50 

0 





| —— C02: calculated 


O CO2: measured 
Propane: calculated 
Propane: measured 














Fig. 12. S-L—V phase equilibria for vitamin K, in CO, or 
propane: comparison of calculated and measured values. 
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Fig. 13. Solubility of vitamin K, in dense gas: comparison of 
experimental and calculated data: (a) CO,; and (b) propane. 


fraction at 40°C and 253 bar. The solubility 
isotherm at 80°C lies on S—L—V three phase line. 

The solubility of vitamin K, in propane (Fig. 
10(b), Table 7) is in the range from 0.025 x 107? 
mole fraction at 40°C and 11.5 bar to 2.99 x 10-3 
mole fraction at 80°C and 65 bar. All isotherms 
present the solubility of a solid vitamin in propane. 
The solubility isotherms run through a concentra- 
tion maxima, which is due to the equality of the 
partial molar volume of the solute in the supercrit- 
ical phase and the molar volume of the pure solid 
substance. The concentration maxima is with in- 
creasing temperature shifted to lower densities 
(Table 7). Similar phenomena were also found for 
some other systems [19,20]. 


Fig. 11 presents the comparison of experimental 
solubilities of vitamin K, in CO, with the literature 
data [1]. It is evident, that the measured solubilities 
are in the same order of magnitude as the literature 
data, but they are somewhat lower. The difference 
is ~ 50-60%. 

It can be concluded that the solubilities of 
vitamins D in supercritical CO, under the process 
conditions investigated are in the range 0.04 x 
10~ ?—1.45 x 10~3 mole fraction, except for vita- 
min K;, where the solubility in CO, is higher and 
is in the range from 0.16 x 10~* to 4.07 x 10-3 
mole fraction. The solubilities of vitamins D in 
propane are approximately up to 10 times higher 
as in CO, and are in the range from 0.15 x 107° 
to 12.40 x 10~* mole fraction. Opposite behaviour 
is observed for vitamin K, where the solubilities in 
propane are lower than in CO, and are in the range 
from 0.03 x 10~* to 2.99 x 10~3 mole fraction. 

According to Tables 5—7 obtained, solubility 
data are in a very wide density range from 164.5 
to 945.5 kg/m? for CO, and 8.9 to 515.9 kg/m? for 
propane. At these conditions propane is in gaseous 
or liquid state and CO, in liquid or SCF state. It 
is known generally, that at constant density the 
solubility increases with increasing temperature. 
However, present data are showing an unusual 
behaviour at lover densities, where solubility starts 
to increase with decreasing temperature and there- 
fore a crossing of isotherms could be observed in 
density plots. In the case of propane such behaviour 
can be observed for vitamins D, and D; at very low 
densities where propane is in gaseous state. In the 
case of CO, 60°C solubility isotherms for K;, D, 
and D; are taking an unusual course, what is 
probably a consequence of solid’s phase transition. 
Similar crossing of isotherms was already observed 
for the system l-amino-4-hydroxy-2 phenoxy-9,10 
antrachinon disperse red (DR60)-CO, by Tuma 
[21], however the reason for crossing the isotherms 
in density plot could not be explained yet. 


3.3. S—L-—V data correlation 


The experimental S-L—V data were correlated 
using Peng-Robinson equation of state and binary 
interaction parameters were obtained from the best 
fit The approach suggested by McHugh and 
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Krukonis [22] for binary systems was used for 
calculation of solid solubilities using binary 
parameters, determined by modelling S—L—V line. 
This approach reduces the amount of experimental 
work and is therefore very interesting. 

The parameters needed in the Peng-Robinson 
equation of state to model S—L—V line are deter- 
mined by solving following equations: 


SUT, P, y2) = f(T, P, Xo) (1) 
AMT, Poy) =f, P, x) (3) 


where f is the fugacity. Superscripts V, L and S 
denote vapor, liquid and solid phase, respectively 
and subscripts 1 and 2 denote light and heavy 
components in the system, respectively. Equation 
2 reduces to: 


(4) 





which presents the solubility of solid component in 
SCF, y, is mole fraction and @, fugacity coefficient 
of component 2 in SCF. ¢§ is the fugacity coeffi- 
cient of component 2 at its vapor pressure and 
equals unity, due to low vapor pressures of solids. 
P§ is vapor pressure and v} mole volume of pure 
solid component, P and 7 are the pressure and 
temperature at equilibrium, R is gas constant. 


Equations 1,3 reduces to: 


V2Pd = X24 (5) 
ViPt =X pr (6) 
The solute properties were estimated by Lydersen 
group contribution methods [23] (Table 8). For the 
calculation of vapor pressure of the solids, the 
approximation, which relates the vapor pressure of 
the solid to that of the supercold liquid, was used 
[24]. The vapor pressure of pure solid can be 
estimated using the fusion properties of the solid 
component: 


AS™S(T.,.) ( T. 
aN Rc i A ie. 7 
In PS =In P} . (2: ) (7) 


where P$ is the solid vapor pressure, P5 is the vapor 
pressure of subcold liquid, AS%'s is the entropy of 


fusion for component 2 at its normal melting point 
T,,, and T is the operating temperature. 

The entropy of fusion for vitamins were calcu- 
lated from fusion enthalpy determined with dif-fer- 
ential scanning calorimetry-DSC  (Thermo- 
analytical system Mettler TA 3000) (Table 8). The 
solid molar volume of the vitamins was determined 
from their densities measured by helium pycnome- 
ter (Table 8). 

The fit of S-L—V lines in the case of CO, and 
propane was performed and binary parameters k, 
and 7,, (Table 9) obtained at optimal fit were used 
for calculations of equilibrium solubilities of vita- 
min in dense gas. The fit of S-L—V lines shows 
good results in the case of vitamin K, (Fig. 12). For 
the vitamins D, and D;, which are substances with 
more complex molecule, the approach used gives 
no good results. The order of the magnitude of 
vitamin K, solubilities is well represented, although 
the data are not quantitatively modelled (Fig. 13). 
The deviation of the model from the data could be 
a consequence of a fact, that binary parameters are 
temperature dependent. Therefore, the solubilities 
can not be described accurately in a wide temper- 
ature and pressure range with single values ob- 
tained from S—L—V fit. 


4. Conclusions 


Phase equilibrium data of binary systems con- 
taining fat-soluble vitamin (D,, D; or K;) and 
dense gas (CO, or propane) were measured in order 
to obtain much needed data for the designing of 
high pressure separation process for isolation of 
vitamins from natural sources. The results have 
shown, that the solubilities of fat soluble vitamins 
D, and D, are approximately up to 10 times higher 
in propane as in CO,. Opposite behaviour was 
observed for vitamin K;, where solubilities in CO, 
are somewhat higher as in propane. In agreement 
with these findings were also the observations of the 
solid—liquid transition of vitamins under the pres- 
sure of gas. The melting point depression of vita- 
mins D, and D; is more pronounced in propane as 
in CO,, whereas the melting point of vitamin K, 
depresses to lower temperatures under the pressure 
of CO, as under the pressure of propane. 
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The results indicate, that higher miscibility of 
two components in binary system consequences in 
higher melting point depression of heavy 
component. 
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Abstract 


A high pressure, variable-volume visual cell was used to perform static measurements of phase equilibria involving 
carbon dioxide and two different fatty acid ethyl esters: eicosapentaenoic acid (EPA) and docosahexaenoic acid 
(DHA) ethyl esters. Six temperatures ranging from 303.15 to 353.15 K were selected. In a second step, our data were 
compared with predictions obtained by two group contribution based models previously published. Both models are 
able to predict with quite high accuracy our data whatever the temperature. © 2001 Elsevier Science B.V. All rights 


reserved. 


Keywords: Fatty acid ethyl esters; Phase equilibria; Solubility; Supercritical CO,; Variable-volume view cell 














1. Introduction 


Separation and fractionation of fish oil related 
compounds have become a subject of active re- 
search because of potential applications in food 
and biochemical industries. In particular EPA 
(eicosapentaenoic acid) and DHA (docosahex- 
aenoic) ethyl esters have been noted for their 
beneficial role in human health and diseases. w-3 
fatty acids are vital for the retina of the human 
eye and the nervous system. They are constituents 


* Corresponding author. Tel.: + 33-383-175081; fax: + 33- 
383-175152. 

E-mail address: jean-noel.jaubert@ensic.inpl-nancy.fr (J.-N. 
Jaubert). 


of our cell membranes, sending signals to 
molecules, hormones and enzymes. -3 fatty acids 
have been shown to reduce the risk of cardiovas- 
cular disease and inflammatory diseases, but they 
also seem to have an influence on mental develop- 
ment and mental health, the duration of preg- 
nancy and sperm quality. 

Due to their thermal sensitivity, however, con- 
ventional separation techniques such as vacuum 
distillation, requiring relatively high temperature, 
are unsuitable. Supercritical fluid (SCF) separa- 
tion processes using carbon dioxide (CO,) as sol- 
vent offers potential advantages for processing 
under milder conditions. Furthermore, CO, is 
nonflammable, inexpensive and non toxic. During 
the actual designing of this kind of supercritical 


0896-8446/01/$ - see front matter © 2001 Elsevier Science B.V. All rights reserved. 
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separation processes, a sufficient knowledge of the 
vapor liquid equilibrium (VLE) data for mixtures 
consisting of fatty acid ethyl esters (FAEE) and 
CO, is indispensable. This is the reason why VLE 
data for binary mixtures of CO, and two FAEE 
are reported in this work. EPA and DHA ethyl 
esters which are both m-3 FAEE were selected 
because of their particular importance mentioned 
above. 

For these two systems, Bharath et al. [1] made, 
> 10 years ago, measurements at three different 
temperatures: 313.15, 323.15 and 333.15 K. How- 
ever their products were not very pure. In their 
paper, the authors indicate that the purity of their 
FAEE was ~ 90%. In this study, we were able to 
work with purer components and in order to test 
the influence of the solute purity on equilibrium 
pressures, measurements were also made at 
313.15, 323.15 and 333.15 K. Moreover, three 
new temperatures were added (303.15, 343.15 and 
353.15 K). 

In a second step, our data were compared with 
predictions obtained by two group contribution 
(GC) based models recently described in the liter- 
ature. Indeed, nowadays, in order to design a 
SCF extractor, it is highly suggested to have at 
one’s disposal purely predictive models which do 
not require the phase equilibrium measurements 
of all the binary systems involved in a complex 
mixture. 

The first model: ‘MHV1-QB approach’, devel- 
oped in 1995 by Coniglio et al. [2], combines the 
Soave-Redlich-K wong [3] equation of state (EOS) 
with the Modified UNIFAC [4] g® model through 
the MHVI [5] mixing rule and with a quadratic 
expression for the mixture b parameter. Because a 
constant /,, value is used, such a model is purely 
predictive. 

The second model, developed by Jaubert et al. 
[6—9], is an extension of Abdoul’s model [10] to 
fish oil related components. Such a model com- 
bines at constant packing fraction the Van Laar 
g® model and the Peng-Robinson [11] EOS. In 
order to have a purely predictive model, the bi- 
nary interaction parameters of the Van Laar func- 
tion are estimated through a group contribution 
method [6-10]. 


In the two previously described models, the 
critical parameters (7. and P.) and the acentric 
factor necessary to solve the EOS are estimated 
through correlations developed by Coniglio et al. 
[2,6] especially for fish oil related components. 


2. Experimental section 
2.1. Material 


Carbon dioxide was supplied by Carboxyque 
Frangaise and was of industrial grade (an average 
purity of 99.998%). EPA and DHA ethyl esters 
were nicely given by Pronova a.s. which is a 
biochemical company owned by Norsk Hydro 
ASA (Pronova Biocare is the World’s largest 
manufacturer of w-3 fatty acids). The purity of 
the two fatty acid ethyl esters (FAEE) were re- 
spectively equal to 95.7% (EPA-EE) and 93.6% 
(DHA-EE). 


2.2. Apparatus and procedure 


Bubble and dew point pressures of mixtures 
involving CO, and a FAEE were obtained using 
the static type apparatus as shown in Fig. 1. Such 
a synthetic method was selected because no sam- 
pling is necessary since a known composition 
mixture is introduced in an optical cell. The phase 
transition resulting from pressure variations are 
studied by direct visualization. By this way, the 
precision of the measurements depends essentially 
on the precision of the weighing. This equipment 
[6,7] consists of an equilibrium cell with a window 
in sapphire. The temperature inside the cell is 
regulated thanks to a thermostated bath. The 
highest working temperature is equal to 150°C. 
The pressure inside the cell may vary between | 
and 700 bars and is measured with an accuracy of 
1 bar thanks to a pressure gauge. The total cell 
volume is recorded with an accuracy of 0.1 ml on 
the graduated piston connected to the cell. Such a 
volume may vary between 1.7 and 12.4 cm?. 

At the beginning of each experiment, the stain- 
less steel screw is removed, i.e. the top of the cell 
is opened and a given amount of FAEE is intro- 
duced into the cell with a syringe. By weighing the 
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syringe before and after the injection, the amount 
of FAEE introduced in the cell can be determined 
very accurately. The cell containing the FAEE 
being opened, the CO, tank is slightly opened in 
order to send CO, inside the cell and to overrun 
air outside. A few seconds later, the cell which 
now contains the FAEE and a small amount of 
CO, is closed thanks to the steel screw. The CO, 
tank is also closed and the cell is thermostated at 
low temperature. Once the temperature and pres- 
sure are uniform inside the cell, the CO, bottle is 
once again opened and a given amount of CO, 
spontaneously enters the cell because of the pres- 
sure difference between the cell and the CO, tank. 
After closing the CO, bottle, the cell is left until 
the pressure and the temperature stabilize. These 
initial temperature and pressure are recorded. The 
initial total cell volume is simply read on the 
graduated piston connected to the cell. Assuming 
that at low temperature the FAEE and the CO, 
are not miscible or more precisely assuming that 
excess volume is equal to zero if mixing takes 
place, it is easy to determine the amount of CO, 
introduced in the cell. To do so, it is necessary to 
know the CO, and the FAEE densities under the 
initial pressure and temperature. The CO, density 
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is know from IUPAC tables [12] whereas FAEE 
density was given by Pronova a.s. Once the global 
cell composition (amount of CO, and FAEE) is 
known, the cell was thermostated at different 
temperatures. For a given temperature, the piston 
was manually moved very slowly in order to 
reduce the total volume and to increase pressure 
until the cell becomes monophasic. In fact, once 
the pressure is stabilized, it is slowly decreased 
until a second phase appears. The medium is thus 
successively compressed and expanded in order to 
define the most narrow range of pressure for the 
phase transition. The stabilization has to be 
reached after each modification of pressure, but 
due to the small cell volume and to the efficiency 
of the stirring, the equilibration does not require 
more than three hours. The pressure gauge was 
thus indicating the saturation pressure of the sys- 
tem, i.e. the pressure at which a given amount of 
matter thermostated at a given temperature be- 
comes monophasic. The experimental accuracy on 
such a measurement was found to be ~1 bar 
since the cine camera connected to the cell allows 
a very precise location of the phase transition. If 
increasing pressure, the liquid phase disappears, a 
dew point pressure is measured. In the opposite 
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Fig. 1. Apparatus used for bubble and dew point measurements: (1) equilibrium cell; (2) window in sapphire; (3) CO, tank; (4) 
removable screwed cap allowing to feed the FAEE; (5) thermostated water outlet; (6) thermostated water inlet; (7) TV connected 
to the cine camera aimed at observing inside the cell; (8) stirring magnet; (9) graduated piston allowing the measurement of the cell 
volume; (10) manual system allowing the pressure to vary by changing the cell volume. 
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Table 1 
Vapor-liquid equilibria for CO,-EPA ethyl ester (e20-5@3) 
system 





Pressure (bar) Mole fraction of CO, 





Liquid phase Vapor phase 





303.15 K 


60.0 : 
64.0 0.9985 


65.0 0.9972 
66.0 0.9997 
70.0 0.9927 


75.0 
79.0 
83.0 
100.0 
109.0 
118.0 


313.15 K 
65.0 
75.0 
77.0 
79.0 
83.0 
95.0 

102.0 

106.0 

120.0 

130.0 

134.0 

139.0 


323.15 K 
80.0 
90.0 
96.0 

104.0 

116.0 

122.0 

129.0 

145.0 

154.0 

163.0 

165.0 


333.15 K 
90.0 
101.0 
110.0 
113.0 
121.0 
131.0 
148.0 
155.0 
169.0 
189.0 


Table 1 (Continued) 





Pressure (bar) Mole fraction of CO, 





Liquid phase Vapor phase 





343.15 K 

102.0 0.7506 _ 
111.0 _ 0.9997 
124.0 0.9985 


126.0 ~ 
141.0 “ 0.9972 


145.0 - 
170.0 0.9927 
173.0 0.9938 
190.0 0.9892 


200.0 o 
212.0 , 0.9833 


353.15 K 
120.0 0.9997 


141.0 _ 
145.0 0.9985 
151.0 0.9972 


162.0 
180.0 
197.0 
210.0 
218.0 
232.0 
233.0 








case (disappearance of the gas phase), a bubble 
point pressure is measured. 


2.3. Experimental results 


A total of 140 saturation pressures (53 bubble 
point pressures and 87 dew point pressures) were 
measured on the systems CO,-EPA ethyl ester 
and CO,-DHA ethyl ester. For each of them, 70 
equilibrium pressures were experimentally deter- 
mined. Concerning more precisely the system 
CO,-EPA ethyl ester, VLE (28 bubble point pres- 
sures and 42 dew point pressures) were measured 
at six different temperatures (T/K = 303.15, 
313.15, 323.15, 333.15, 343.15 and 353.15). Con- 
cerning the system CO,-DHA ethyl ester, VLE 
(25 bubble point pressures and 45 dew point 
pressures) were measured at five different temper- 
atures (T/K = 313.15, 323.15, 333.15, 343.15 and 
353.15). The experimentally measured data are 
given in Tables 1 and 2 and are visible in Figs, 2 
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Fig. 2. Experimental and calculated vapor—liquid equilibria for the carbon dioxide-EPA ethyl ester system. M and @, bubble and 


dew point pressures previously measured by Bharath et al. [1]; (J and O, bubble and dew point pressures measured in this study; 


——., predicted phase envelope using Jaubert et al.’s model [6—9] (extension of Abdoul’s model [10] to fish oil related components); 
, predicted phase envelope using the MHV1-QB approach (Coniglio et al. [2]). (a) 7/K = 303.15; (b) T/K = 313.15; (c) 


T/K = 323.15; (d) T/K = 333.15; (e) T/K = 343.15; (f) T/K = 353.15. 
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Fig. 3. Experimental and calculated vapor—liquid equilibria for the carbon dioxide DHA ethyl ester system. M@ and @, bubble and 
dew point pressures previously measured by Bharath et al. [1]; (1 and ©, bubble and dew point pressures measured in this study; 
——., predicted phase envelope using Jaubert et al.’s model [6—9] (extension of Abdoul’s model [10] to fish oil related components); 
—~—-—- ~~, predicted phase envelope using the MHV1-QB approach (Coniglio et al. [2]). (a) T/K = 313.15; (b) T/K = 323.15; (¢) 


T/K = 333.15; (d) T/K = 343.15; (e) T/K = 353.15. 
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Table 2 
Vapor—liquid equilibria for CO,-DHA ethyi ester (e22-63) 
system 











Pressure (bar) Mole fraction of CO, 





Liquid phase Vapor phase 














313.15 K 
47.0 0.6477 
69.0 0.7717 
80.0 
81.0 
84.0 
92.0 
93.0 
110.0 
116.0 
116.0 
117.0 
136.0 
146.0 
163.0 
175.0 


323.15 K 
55.0 
83.0 
97.0 
99.0 

100.0 

114.0 

120.0 

135.0 

136.0 

144.0 

147.0 

158.0 

164.0 

172.0 

184.0 

200.0 


333.15 K 
63.0 
102.0 
109.0 
115.0 
120.0 
130.0 
143.0 
158.0 
166.0 
169.0 
181.0 
190.0 
197.0 
210.0 
225.0 


Table 2 (Continued) 





Pressure (bar) Mole fraction of CO, 





Liquid phase Vapor phase 





343.15 K 
75.0 0.6477 
121.0 0.7717 
125.0 - 
132.0 
138.0 
139.0 
162.0 
175.0 
188.0 
190.0 
211.0 
221.0 


353.15 K 
96.0 
134.0 
144.0 
148.0 
154.0 
189.0 
196.0 
207.0 
211.0 
231.0 
242.0 





and 3. In these two figures are also given the 
data points previously measured by Bharath et 
al. [1]. 

From these two figures, it is quite difficult to 
compare our data with those previously deter- 
mined by Bharath et al. [1]. This means, it is 
difficult to appreciate the benefices obtained by 
using purer solutes. At the lowest temperatures, 
i.e. at 313.15 and 323.15 K, our measured bub- 
ble point pressures are higher than those of 
Bharath et al. except in the vicinity of the criti- 
cal point where the reverse situation occurs. At 
the highest temperature, i.e. at 333.15 K, the 
both set of data are in good agreement. For the 
CO,-DHA ethyl ester system, our estimated crit- 
ical point seems however lower than the one 
estimated from Bharath’s data. 
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3. Prediction of the data 


As stated in the introduction, two GC based 
models were used in order to predict our data. 
These two models are purely predictive in the 
sense that all the parameters required in the equa- 
tion of state (EOS) are estimated with the unique 
knowledge of solute (EPA EE or DHA EE) 
molecular structure. This class of models has the 
advantage to make time and money saved by 
avoiding the use of any measurement and are 
therefore particularly appreciated by private com- 
panies. As an example, it is not necessary to 
perform VLE measurements for all the binary 
systems included in a complex mixture in order to 
estimate its phase behavior. 


3.1. Coniglio et al.’s model (MHV1-QB approach) 


The first model we used was developed in 1995 
by Coniglio and coworkers [2]. In this approach, 
all the calculations are performed by means of the 
Soave-Redlich-Kwong (SRK) equation of state 





(1) 


where a and b are the mixture parameters. 
A quadratic expression was selected for the 
mixture b parameter according to Equation 2. 


a (O55 + bj) 


b=) 2) ux 5) 


i=1j=1 


(i —/,) (2) 


where 5, and 5, are the pure component b 
parameters calculated by: 


ds oe 
and og, «v2 (3) 


RT. 
b..= a 
li Q, P 3 


Sj 


Table 3 


Furthermore, in order to obtain a purely predic- 
tive model, Coniglio et al. chose a constant value 
for the binary interaction parameter /,, (/;, = 0.3). 

The mixture a parameter is estimated from the 
first order Modified Huron-Vidal mixing rule 
(MHV1) proposed by Michelsen [5]: 


a ee a sie (2 
——_ = Zz _ eee x1 a, 
eS MI aie ibe 3h ie 2 
with g,= —0.593 (4) 


In equation 4, Modified UNIFAC [4] was selected 
as the g® model (new values of the CO,-solute 
group interaction parameters were defined by 
Coniglio et al.). a;; is the pure component attrac- 
tive parameter calculated by: 
ee be a Pe 4 

P. [ATP with T, = T. (5) 


1 


where f(T, ) is the expression defined by Mathias 
and Copeman [13]. 


qij = ().4286 


AT.) = 


eee ee ee 


l+C(1-/T,) if T,>1 


The values of the constants C,, C3, C; and those 
of the critical parameters (7., P.) estimated 
through the Ambrose [14,15] group contribution 
methods (GCM) are given in Table 3 for the EPA 
and DHA ethyl esters. The normal boiling tem- 
perature given in Table 3 which is required for the 
Ambrose GCM was estimated using the AMP 
equation [16]. Such an equation allows the calcu- 
lation of the vapor pressure of pure components. 
Its parameters were estimated by the correlations 
developed by Coniglio et al. [2]. By the end, by 


Estimated pure components physical properties and C; parameters of the Mathias and Copeman [13] expression 





Compound T,/K T./K P.,/bar 


@M C; Cc. C; 





792.66 12.06 
803.10 11.28 


640.81 
655.63 


EPA ethyl ester 
DHA ethyl ester 


0.1207 
— 0.1454 


— 0.3566 
— 0.1943 


1.9562 
2.0243 


1.0102 
1.0797 
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knowing the critical temperature and the AMP 
equation parameters, the acentric factor (not 
needed in MHVI1-QB approach) was simply esti- 
mated by: 


o= — log, )P%*"( 2 — 0.7) — l (7) 


The prediction of our data with this first model 
are visible in Figs. 2 and 3 (dashed lines). A 
bubble point pressure and a dew point pressure 
algorithms were used to build the curves. From 
these figures, it is possible to conclude that such a 
model is particularly suitable at low pressures but 
overestimates the pressures in the critical region. 
Also, its predictive power is better at high 
temperatures. 


3.2. Jaubert et al.’ s model (extension of Abdoul’s 
model to fish oil related components) 

The second model we used was developed in 
1999 by Jaubert and coworkers [6—9]. In this 
approach, all the calculations are performed by 
means of the Peng-Robinson (PR) equation of 
state [11]: 

RT a 
a ae ©) 
)+ b(v—b) 
To describe the behavior of pure components, the 
following Soave [3] type function, modified by 
Rauzy [17] was used: 


R27? T. \0-445075\ 72 
, = 0.457236 <i : i? foe pee 
aremmmereinlinte)  ) 


. (9) 





with: 


m,; = 
6.812553[,/ 1.127539 + 0.517252, — 0.003737? — 1] 
(10) 








The covolume of each pure component is classi- 
cally calculated by: 


b, = 0.0778 0-2: (11) 


Cc. 
I 


The pure component physical properties re- 
quested to use such an EOS (T7., P., @) are 
estimated exactly as in the previous model pro- 
posed by Coniglio et al. (see Table 3). 


153 


In Jaubert et al.’ s approach, the mixing rules 
for estimating a and b mixture parameters are 
defined by introducing at constant packing frac- 
tion an excess Gibbs energy function (g*). From 
classical thermodynamics, it is well known that 
the fugacity coefficient gy; of a component / in a 
homogeneous phase is given by: 


AA res RT 
In g;= sel baad —Inz (12) 
On; Eis Ni zi 


where z is the compressibility factor of the phase. 

By defining the packing fraction of a phase by 
yn =b/v where b and v are respectively the covol- 
ume and the molar volume of the phase, one 


obtains: 





A(T, , n) ue 
ee 





y » ae 
n| Sol ac em (13) 
0 ] 


In order to compute the fugacity coefficients 
through equation 12, it is thus necessary to be 
able to estimate the compressibility factor z(7, y, 
n) of a phase. According to the reference at 
constant packing fraction, Péneloux et al. [18] 
have shown that z(7, 7, n) could be expressed by: 


Pp 
2(T,9,x)= ¥ x,2¢T, 9) + 277, 9, x) (14) 


i=1 


with: 





E 
n (= (T, 7, ») (15) 
1.x 


-E tices ea 
a 1 O= Re On 

z, is the compressibility factor of pure component 
i calculated by the PR-EOS. Moreover, in equa- 
tion 15, A may be expressed as explained by 
Péneloux et al. [18] as the product of two func- 
tions. The first one is a classical excess Gibbs 
energy function (Van Laar, NRTL, UNIQUAC) 
depending only on temperature and composition, 
whereas the second one only depends on the 
packing fraction. This means: 


A®(T, n, x) =g"(T, x) x F() (16) 


According to the previous work of Abdoul et al. 
[10], the following function was selected for F(7): 
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Table 4 
Deviations observed by predicting our data with Jaubert’s approach 








System: CO,+ T/K N,* AP% 


b 
Np. bubble 





EPA ethyl ester 303.15 12 12.00 
313.45 12 10.69 
323.15 12 8.30 
333.15 11 10.08 
343.15 12 10.75 
353.15 11 10.01 


Global 303-353 70 10.31 


DHA ethyl ester 313.15 15 9.78 
323.15 16 7.85 
333.15 16 7.03 
343.15 12 5.60 
353.15 1] 4.86 


Global 313-353 70 7.22 25 


Both systems 303-353 140 8.77 53 








* N,, is the number of experimentally determined saturation pressures (bubble + dew). 
> Np.bubbie 1S the number of experimentally determined bubble point pressures. 
© N,.dew 18 the number of experimentally determined dew point pressures. 


pared to our data, the deviations observed be- 
tween calculated and experimental saturation 
(bubble and dew) pressures are given in Table 4. 
The global average deviation observed for the two 
systems investigated (140 measured saturation 
pressures) is equal to 8.77%. Such a deviation 
et with lower than 10% is completely acceptable for an 
2 industrial use and shows the interest of purely 
predictive thermodynamic models. 





| wi+n(l +/2) on 


F(y)=— 
” 2/2 1+n(01—./2) 


and a Van Laar like g™ function was selected: 


7 a 
Z Zz X;X; 5; Ey(T) 


binkin4 





g"(T, x) 


(18) 


E,(T) 22 E\(T) 
E,(T) =0.0 


In equation 18, the covolume b, of each pure 
component i is calculated through equation 11 
whereas the binary interaction parameter £,(7) 
was computed by the GCM initially developed by 
Abdoul et al. [10] for hydrocarbons and extended 
to oxygenated components by Jaubert et al. [6—9]. 

The prediction of our data with this second 
model are visible in Figs. 2 and 3 (continuous 
lines). From these figures, it is possible to con- 
clude that such a model is particularly suitable in 
the critical region. Moreover, the dew point pres- 
sures are predicted very accurately. Because this 
second model leads to better prediction as com- 


4. Conclusion 


A high pressure, variable-volume view cell was 
used to perform static measurements of equilibria 
for two binary systems involving CO, and a 
FAEE. The two selected solutes were EPA and 
DHA ethyl esters because of their importance in 
food, cosmetics and biochemical industries. For 
each system, 70 equilibrium pressures (dew or 
bubble) were measured. Six temperatures ranging 
from 303.15 to 353.15 K were selected. 

Our data were compared with predictions ob- 
tained by using two GC based models recently 
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described in the literature. The first one developed 
by Coniglio et al. in 1995 leads to very good 
results at small pressures but overestimates the 
bubble point pressures in the vicinity of the criti- 
cal point and all the dew point pressures. The 
second one developed by Jaubert et al. (1999) was 
able to correlate the 53 bubble point pressures 
with an average overall deviation of 11.09% and 
the 87 dew point pressures with an average overall 
deviation of 7.35%. This means that the 140 satu- 
ration pressures were predicted with an average 
deviation of 8.77% what is really very satisfactory 
for industrial applications. 

The better performance of Jaubert et al.’s 
model in the vicinity of the critical point and for 
the dew point prediction as compared to 
Coniglio’s model which has better performance 
under low pressure, can be explained by some 
intrinsic features of the two models: (1) MHV1- 
QB is a mixing rule where zero pressure has been 
chosen as reference. The constant packing frac- 
tion taken as reference in the Abdoul’s model is 
more or less equivalent to a reference pressure 
which would be variable depending on the mix- 


ture composition and temperature; (2) model GC 
parameters were regressed on binary VLE data by 
using a (T, P) flash algorithm for the approach of 
Jaubert et al., whereas a bubble point algorithm 
was used in the approach developed by Coniglio 


et al. 
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Abstract 


Supercritical fluid extraction posses several advantages Over traditional liquid—solvent-based extraction methods 
including improved selectivity, expeditiousness, automation and environmental safety. Fluid phase equilibrium data 
is essential for study the viability of supercritical extraction and the design of the extraction columns often used for 
this type of separations. During the past few years it has been demonstrated that phenolic compounds, such catechins 
and procyanidins from seeds, contribute to both the flavour properties and the pharmacological effects of wine. The 
extraction of catechin using supercritical CO, is not possible, it only extracts traces of product, probably due to its 
polar nature. So, it is necessary the use of a cosolvent. Ethanol is appropriate because it is a polar solvent permitted 
in food industry. This work uses different mixtures ethanol—CO,, ranging 5—30% the ethanol contained, to optimise 
the amount of catechin extracted. The solubility of catechin in their supercritical solvent was measured at 313 K and 
pressures ranging 80-120 bar. In order to correlate the obtained results, two equations of state (EOS) were 
successfully used (Peng—Robinson (PR) and Soave—Redlich- Kwong (SRK)). Generally both equations show similar 
deviations, but PR shows the better correlation with experimental data. The solubility of catechin shows an increase 
with the value of percentage of ethanol. The solubility rises with pressure. The extraction pressure will be limited by 
security cautions and, just such ethanol percentage, by operational costs. © 2001 Elsevier Science B.V. All rights 
reserved. 


Keywords: Supercritical solubility; Catechin; Carbon dioxide; Ethanol; Modelling 














highly valuable natural products. Supercritical 
fluid extraction posses several advantages over 
traditional liquid—solvent-based extraction meth- 
ods including improved selectivity, expeditious- 
ness, automation and environmental safety [1—4]. 
cle tdee The avoidance of organic solvents is a major goal 

* Corresponding author. Tel.: + 34-96-3864316; fax: + 34- in the isolation of natural products, which may be 
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1. Introduction 


There is an urgent need for rapid and clean 
methods for the extraction and determination of 
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During the past few years it has been demon- 
strated that phenolic compounds contribute to 
both the flavour properties and the pharmacologi- 
cal effects of wine. Phenolic compounds such 
catechins and procyanidins from seeds have been 
shown to affect the bitterness and astringency of 
wines [5—7]. These phenols may act against in 
vitro oxidation of low-density lipoprotein [8]. It 
also has been suggested to have antiulcer [9], 
anticarcinogenic [10], antimutagenic [11], and an- 
tiviral [12] activities. The high antioxidant power 
of phenols is generally believed to account for 
these activities [13,14]. They also exhibit various 
physiological activities including antiinflamma- 
tory, antiallergenic, antihypertensive and an- 
tiarthritic activities [15]. 

Common methods for the isolation of phenolic 
compounds from grape seed use organic solvents 
such methanol [16-18], ethanol [18,19] and ace- 
tone [20]. Extraction recoveries improve when 
mixtures of methanol—water and acetone—water 
are employed [18,20]. Various extraction tempera- 
tures have been used with times ranging from a 
few minutes to several hours. The most general 
methodology involves the use of aqueous 
methanol for 16—24 h at room temperature and 
atmospheric pressure. 

Supercritical fluid extraction (SFC) affords big 
advantages over more conventional extraction 
techniques. The absence of both light and air 
during the extraction process can reduce the inci- 
dence of degradation reactions that can more 
easily occur with other extraction techniques. For 
example, Tipsrisukond et al. [21] have reported 
that supercritical fluid extracts showed higher an- 
tioxidant power than extracts obtained by classi- 
cal methods. Furthermore, it is not necessary to 
concentrate the resulting extract because only a 
small volume of organic solvent is normally used. 
SFE has been previously applied to grape seeds 
for the removal of oils [22]. Pure CO, without an 
organic modifier was used in this extraction. Phe- 
nolic compounds were apparently not removed 
under these conditions. The use of small quanti- 
ties of organic solvent with supercritical CO, facil- 
itate the phenols extraction from the seed [23-26]. 

High-pressure vapor—liquid equilibrium data 
for CO, + catechin system are required for ratio- 


nal equipment design. There are some tests of 
supercritical extraction of phenolics compounds 
including catechin in literature [23—27]. These au- 
thors study the extraction of phenolics from inert 
supports or the fractionation of extracts from 
white grapes without obtaining any conclusion 
about the viability of supercritical extraction of 
these phenolics. The first step in order to obtain 
this information is the study of the solubility of 
these compounds in supercritical CO,. Literature 
data show the extraction of catechin using CO, 
supercritical was not possible [23-27]. It only 
extracts traces of product. That is probably due to 
the catechin polar nature. So, it is necessary the 
use of a cosolvent. Ethanol is appropriate because 
it is a polar solvent permitted in food industry. 

This work has as objective the determination of 
solubilities of catechin in a mixture ethanol + 
CO,. This is a previous step for the study of the 
viability of supercritical extraction of phenolics 
(mainly catechin, epicatechin and quercetin) from 
skin grape in order to increase its presence in 
wines or to pharmacological use. Also, this infor- 
mation is necessary to optimise operating condi- 
tions of the extraction. 


1.1. Models 


In order to correlate the obtained results, the 
most usual equations of state will be tried (Peng— 
Robinson [28] Eq. (1) and Soave—Redlich— 
Kwong (Soave [29] Eq. (2)). 


RT a(T) 
pe b vy*+2bv—b* 


RT a(T) 
v—b viv+b) 


Where v is molar volume, a (it is a measure of 
intermolecular forces) and bd (it is concerned with 
size of rigid spheres) are parameters of equations, 
P is pressure, TJ is temperature and R is the 
universal constant of gases. 

These models were applied by using the critical 
properties and acentric factor of pure compo- 
nents, and the following conventional mixing 
rules: 


ay = / 4,4; (1 — ky), (3) 
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b,+ 5, 
ie Be 6 oe 1 i;)- (4) 

Where a,, a,b, b; are the parameters of pure 
components (function of critical properties and 
acentric factor) and k,,, ,; are binary adjustable 
interaction parameters. Another mixing rules were 


tried without improving results. 


by = 


2. Experimental section 
2.1. Chemicals 


Catechin (98 mass%, GC grade) was purchased 
from Aldrich Chemical Co. Inc. Ethanol absolute 
(99.8 mass %, GC grade) supplied by Prolabo 
S.A. was used as modifier and solvent to collect 
the extract. The reagents were used without fur- 
ther purification after chromatography failed to 
show any significant impurities. High-purity CO, 
(more than 99.9 vol% purity, SFC grade) was 
used as received. Silica (more than 99.8 mass”) 
extra pure (from Sigma Chemical Co. Inc.) was 
utilised as support of products. The properties of 
the pure components are listed in Table 1. 


2.2. Equipment and procedure 


The extractions were realised by an SFX 3560 
Extractor with two Model 260D syringe pumps 
manufactured by ISCO (Lincoln, NE, USA). 
With the proper plumbing, two pump system can 
deliver modified supercritical fluid or a continuos 
flow of supercritical fluid. The cylinder capacity of 
the pump is 266 ml and the maximum pressure 
510 atm. The temperature can ranging from 313.2 
to 423.2 K, and supercritical fluid flow rate be- 


Table 1 
Critical constants and acentric factor of pure components 








Substance T.. (K) P.. (bar) 





Co, 304.28 73.78 
Ethanol 513.928 61.48 
Catechin 945.0° 46.8> 








* Daubert and Danner ([32]). 
> Estimated by the Constantinou and Gani method, ([33]). 


© Estimated by the Han and Peng method, ([34]). 


tween 0.5 and 5 ml/min. The SFX 3560 sample 
reel holds up to 24 sample cartridges and collec- 
tion vials. The software supplied with the extrac- 
tor uses a dynamic method programming scheme 
which allows entry changes to methods while a 
previously programmed method is running. The 
schematic diagram of the apparatus is shown in a 
previous work [30]. 

The extractor measured the amount of super- 
critical CO, and modifier ethanol used during the 
extraction. The flow rate of supercritical fluid 
used at all experiences was 2 ml/min. This is the 
suitable flow to achieve the equilibrium with satu- 
ration of supercritical phase. Ethanol was used to 
trap the catechin extracted. After that, the extract 
was transferred to 10 ml calibrated flasks and 
diluted with ethanol. In some cases, a higher 
dilution than this could be necessary to achieve a 
concentration level adequate for UV spectrometry 
determination. 

For sample preparation, first 0.5 g of catechin 
was dissolved in 25 ml of ethanol. Second, 2 ml of 
solution was added of diatomaceous earth (1 g) 
previously weighed in the sample cartridge. Fi- 
nally, the spiked solid was allowed to heat (60°C) 
for at least 2 h in order to evaporate ethanol and 
precipitate catechin over the surface of solid parti- 
cles. This procedure is necessary in order to avoid 
preferential ways during extraction and the for- 
mation of masses of particles. 

A 8453 Hewlett-Packard diode array UV—VIS 
spectrophotometer was used for the determina- 
tions of catechin in the extract, after calibration 
with gravimetrically prepared standard solutions. 

Solubility is measured as the slope of catechin 
moles extracted vs. solvent moles (CO, + ethanol) 
at different extraction times (five values at least). 


3. Results and discussion 


The extraction of catechin using CO, supercriti- 
cal was tested. At extreme conditions of pressure 
and temperature (400 bar and 100°C) it only 
extracts traces of product (< 10~° mole fraction) 
in accordance with Tena et al. [23]. That can be 
due to the polar nature of catechin. So, it is 
necessary to use a cosolvent. Ethanol is appropri- 
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Table 2 

Solubilities (10* y,) of catechin (2) expressed as mole fraction, 
in supercritical phase (CO, (3)+ethanol (1)) at 313 K and 
different percentages of modifier 





P (bar) “% of modifier (ethanol (1)) 





25 30 





9.68 15.09 
32.06 





ate because it is a polar solvent and it is permitted 
in food industry for its harmlessness. Also, it is 
easily eliminated of the extract by evaporation at 
room temperature. 

The solubilities of catechin in supercritical 
phase were measured at 313 K, pressures ranging 
80—120 bar and different mixtures ethanol—CO,. 
These results are shown in Table 2. 

This work uses different mixtures ethanol—CO,, 
ranging 5—30% of ethanol, in order to optimise 
the amount of catechin extracted with the mini- 
mal consumption of cosolvent. So, the extract is 
cleaner and the cost of solvent and separation is 
smaller. 

Solubilities were measured at 313 K. This is the 
minimal temperature allowed for the equipment. 
This temperature has been selected in order to 
avoid the opposed effect of temperature on solu- 
bility [1]. 

The inferior limit of 80 bar was chosen because 
it was near of CO, critical pressure, so the be- 
haviour of system at conditions near the critical of 
system was studied. It was not possible to mea- 
sure the solubility of catechin at pressures up of 
110 bar with this equipment because, at these 


Table 3 


conditions, the entire product contained in the 
sample cartridge was extracted. We can measure it 
just to 130 bar only using 5% of ethanol. 

In order to correlate the obtained results, two 
equations of state were successfully used: Peng— 
Robinson (PR EOS) [28] and Soave—Redlich— 
Kwong (SRK EOS) [29]. The parameters of 
equations of state were determined using the 
solver function of Microsoft Excel and PE pro- 
gram [31]. These parameters and the average rela- 
tive deviation (ARD) are shown in Table 3. The 
parameters estimation was performed considering 
all the experimental data by minimising the fol- 
lowing objective function: 


ABS exp ‘ca o) 
= (y p Veal (5) 
exp 


OF = 





Where y is the solubility of solute (catechin) in 
supercritical phase. 

Generally both equations show similar devia- 
tions at different conditions of pressure and tem- 
perature, but PR EOS [28] shows the better 
correlation with experimental points. 

The evolution of the solubility of catechin vs. 
the modifier fraction shows an increase with the 
amount of ethanol employed. So the optimal con- 
centration of cosolvent will be fixed by separation 
costs. This behaviour was shown in Fig. |. In a 
previous work [30], the optimal concentration of 
ethanol for the resveratrol supercritical extraction 
is founded. This concentration was 7.5%. This 
different behaviour (Fig. 2) for these two solutes 
can facilitate their obtaining by fractionation. 

The solubility, as it was shown in Fig. 1, rises 
with pressure. The extraction will be carried out 
at elevated pressure. The extraction pressure will 
be limited by operational costs and _ security 
cautions. 


Parameters and deviations between experimental and calculated values of system ethanol (1)+catechin (2)+CO, (3) at 313.15 K* 








Model 


12 








PR EOS 0.031 
SRK EOS 0.052 


— 0.0172 
— 0.183 








* ARD = 1/N x [Vex —Veatl/Yex: 
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Fig. 1. Solubility (10* y,) of catechin (2) in supercritical CO, 
(3) at 313 K at different % of modifier ethanol (1) of: (MH) 5%; 
(A) 10%; (@) 15%; () 20%; (CJ) 25%; (A) 30%. (—) 
Modelling by PR EOS ([27]). 


As conclusions, the present study reveals the 
influence of process conditions on the solubility of 
catechin in supercritical CO,. It was necessary the 
use of a cosolvent and ethanol was appropriate. 
The solubility increases with the amount of etha- 
nol. It is possible the fractionation of the catechin 
and resveratrol extracts by supercritical extraction 
by increasing the ethanol concentration. The solu- 
bility rises with pressure, so the extraction will be 
carried out at elevated pressure. Finally, PR EOS 
[28] was successfully applied to the ternary system 
catechin—ethanol—CO,. 








% ethanol 


Fig. 2. Solubility (10* y) at different % of modifier ethanol (1) 
of: (@) catechin (2) in supercritical CO, (3) and (1) resvera- 
trol in supercritical CO, (3) at 313 K and 90 bar. 
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Abstract 


The equilibrium solubility of uranyl nitrate in supercritical fluid carbon dioxide (SF-CO,) in the absence and 
presence of different modifiers was investigated. In the presence of methanol and HNO,, uranyl nitrate showed a 
considerable solubility in SF-CO,. SF-CO, modified with different co-extractants such as tributylphosphate (TBP), 
methylisobutyl ketone (MIBK), acetylacetone (AA) and methanol (MeOH) at 60°C and 250 atm was used for the 
extraction of uranyl nitrate from cellulose-based filter papers. The extraction efficiency of the system for uranyl ion 
was found to decrease in the order TBP > MIBK > AA > MeOH. The use of eight potential chelating agents 
bis(2-ethylhexyl)hydrogen phosphate (HDEPH), tri-n-octylphosphine oxide (TOPO), dicyclohexyl-18-crown-6 
(DC18C6), dibenzoylmethane (DBM), 8-hydroxyquinoline (HOX), diphenylaminesulfonic acid (DPASA) 1,4-bis-[4- 
methyl-5-phenyl-2-oxazolyl]benzene (ODMPOPOP) and TBP for uranyl ion extraction from filter paper samples by 
methanol modified SF-CO, was investigated. The organophosphorous reagents TOPO and HDEPH showed the 
highest extraction efficiency in the series. The synergistic extraction of uranyl nitrate with HDEPH and TOPO in the 
presence of DC18C6 was also studied. © 2001 Elsevier Science B.V. All rights reserved. 


Keywords: U(VI); SFE; Modifier; Chelating agent 





temperature—pressure dependence of solvent 
strength are among the attractive properties which 
make supercritical fluids excellent candidates for 
the extraction and recovery of organic compounds 
from solid matrials [4—6]. The most widely used 
supercritical fluid is carbon dioxide because it is 
nontoxic, nonflammable and relatively inexpen- 
sive and possesses reasonable critical properties as 
well as high solvent power for a wide range of 
nonpolar and intermediately polar organic 


1. Introduction 


In the past two decades increasing attention has 
been directed on the supercritical fluid extraction 
(SFE) as an interesting alternative to conventional 
extraction methods [1-3]. The high diffusivity and 
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In recent years, the supercritical fluids modified 
by addition of complexing agents have been used 
in the development of analytical methods for the 
extraction of metal ions from various solid and 
liquid matrices [7]. The most critical aspect of 
such chelation-SFE of metal ions is the selection 
of suitable chelating agents. Important character- 
istics include high stability and fast kinetics of the 
metal complexes, high solubilities of the chelating 
agents and their metal complexes in supercritical 
CO, and complexing selectivity of the chelating 
agents [8]. Wai and co-workers have reported the 
successful use of fluorinated dithiocarbamates and 
B-diketones as effective chelating agents for ex- 
tracting several transition metal ions from solid 
and aqueous samples by SFE [8-14]. They have 
also described the SFE of uranium and thorium 
from nitric acid solutions with some oragnophos- 
phorous reagents [15]. 

Extraction techniques for the separation and 
recovery of uranium from nuclear wastes are of 
critical importance. Liquid—liquid extraction of 
uranium with organic solution of tri-n-octylphos- 
phine oxide [16,17], tri-7-butyl phosphate [18,19], 
crown ethers [20,21], tri-n-octylamine [22,23] or 
bis(2-ethylhexyl)sulfoxide [24] has attracted con- 
siderable attention. However, these classical ex- 
traction methods are usually time-consuming and 
labor-intensive and require large volumes of high- 
purity solvents. In order to reduce the consump- 
tion of and exposure to organic solvent, disposal 
costs, and extraction time, the solid phase extrac- 
tion of ultra trace uranium(VI) using octadecyl 
silica membrane disks modified by TOPO has also 
been reported [25]. 

We have been recently involved in solubility 
studies of dihydroxybenzene isomers [26] and 
some analytically important anthraquinone 
[27,28] and hydroxyxanthone derivatives [29] in 
supercritical carbon dioxide. In this paper, due to 
the increasing interest in the elimination of trace 
amounts of radioactive chemicals from biological 
matrices, we were interested to study the super- 
critical fluid extraction of uranium from cellulose- 
based filter papers, as a model for biological 
samples such as plants and animal tissues, in the 
presence of several different co-extractants. 


2. Experimental 
2.1. Reagents and materials 


The chelating agents tri-n-octylphosphine oxide 
(TOPO), tri-n-butylphosphate (TBP), bis(2-ethyl- 
hexyl)hydrogen phosphate (HDEPH), 1,4-bis-[4- 
methyl-5-phenyl-zoxazolyl]benzene (DMPOPOP), 
8-hydroxyquinoline (HOX), acetylacetone (AA), 
diphenylaminesulfonic acid (DPASA), dibenzoy!l- 
methane (DBM) and dicyclohexyl 18-crown-6 
(DC18C6) were purchased from Merck chemical 
company and used without any further purifica- 
tion except for vacuum drying. HPLC-grade 
methanol (Aldrich) was used as received. All acids 
used were of analytical reagent grade from Merck 
chemical company. All other reagents and sol- 
vents used in this study were of the highest purity 
available from Merck and Aldrich chemical com- 
panies. Pure carbon dioxide (Sabalan, Tehran, 
99.99%) was used for all SFEs. 

A stock solution of uranyl ion (2000 pg ml~') 
was prepared by dissolving an appropriate 
amount of UO,(NO;),°6H,O in doubly distilled 
deionized water in a 100 ml flask. Samples for 
SFE experiments were prepared by spiking 50 ul 
UO3* solution on prewashed filter papers (What- 
man No. 1, 3.0 cm, Cat. No. 1001 030, washed 
with HNO, and rinsed with water). The spiked 
filter papers were allowed to air dry at room 
temperature. 


2.2. Apparatus and procedure 


2.2.1. Solubility studies 

A Suprex MPS/225 integrated SFE—SFC sys- 
tem equipped with a modified static system [26— 
29] for the solubility determination in SFE mode 
was used. Solubility measurements were accom- 
plished with a | ml extraction vessel in the pres- 
sures of 150 and 250 atm at 60°C. The solid 
solutes (100 mg) were mixed well with some 1 g of 
glass beads and packed into the extraction vessel. 
Sintered stainless steel filters (5 um) were used to 
prevent any carry over of the solutes. The static 
equilibrium time was determined by varying the 
extraction time until a constant solubility value 
was obtained; 30 min was thus used as the equi- 
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librium time. Details of the solubility measure- 
ments are given elsewhere [26—28]. 


2.2.2. SFE studies 

In SFE experiments, a home-built stainless steel 
restrictor with appropriate flow rate (0.8 ml 
min~! at 250 atm and 60°C) was used to main- 
tain the extraction pressure. In order to prevent 
sample plugging, the restrict point was warmed 
electrically. The extract was sampled in a collec- 
tion vessel (trapped in 5 ml methanol). Determi- 
nation of uranium  was_ carried out 
spectrophotometrically using a Model 2100 Shi- 
madzu UV is spectrophotometer. 

A two-step extraction process including 15 min 
static and 15 min dynamic extraction was used. 
After introducing the desired amount of ligand 
into the prespiked filter paper, it was placed into 


Table | 
Solubility of naphthalene in supercritical carbon dioxide at 
35°C 





P(bar) 10°» 








This work [30] [31] 





107 11.6 
138 15.0 
168 16.3 
204 17.5 








Table 2 
Solubility of uranyl ion in pure and modified supercritical 


carbon dioxidie 





Modifier Solubility (g 1~') 








60°C, 150 atm 60°C, 250 atm 





None 0.015 + 0.002 

50 pl MeOH 2.67 + 0.83 

50 wl MeOH +50 pl 0.111 + 0.078 
H,O 

50 wl MeOH +50 ul 
of 1 M HNO, 

50 wl MeOH +50 ul 
of 1M 
CH;,COOH 

50 pl MeOH +50 pl 
of IM HCl 


0.018 + 0.003 
5.26 + 0.24 
0.175 + 0.012 
43.6 + 2.3 61.1 + 3.3 


43.94+5.7 


0.093 + 0.020 








the extraction cell and was then spiked with ap- 
propriate amount of modifier (solvent modifier or 
matrix modifier). Then, the extraction fluid was 
delivered to the extraction vessel (load position). 
After 15 min static extraction, the outlet valve 
was opened automatically and the sample was 
extracted dynamically with SF-CO, for 15 min at 
the same pressure and temperature. In the dy- 
namic extraction step, a supercritical fluid flow 
rate of 0.8 ml min~ ! was used, in order to prevent 
the restrictor plugging. The depressurized fluid 
was bubbled through a collection vial containing 
5 ml of methanol (collect solution). After the 
dynamic extraction step was completed, the filter 
paper was removed from the cell using a pincer 
and placed in a volumetric flask containing 5 ml 
methanol and the flask was shaken vigorously to 
extract the UO3* ions from the filter paper into 
methanol (wash solution). The filter paper was 
finally placed on the neck of the volumetric flask 
and washed with another 2 ml portion of 
methanol. To the ‘collect’? and ‘wash’ solutions 
was then added | ml of dibenzoylmethane solu- 
tion (2% w/v in 50% pyridine—methanol) and 
diluted to the mark with methanol. The uranyl 
concentration was finally determined by ab- 
sorbance measurement of the ‘collect’ and ‘wash’ 
solutions at 405 nm against a reagent blank, using 
an appropriate calibration curve (linear range of 
0.2—6.0 ng ml~' with r= 0.9999). 


3. Results and discussion 
3.1. Solubilily of uranyl nitrate in SF-CO, 


The reliability of the apparatus was preliminary 
checked by measuring the solubilities (mole frac- 
tion, y) of naphthalene in supercritical CO, at 
35°C and the results are listed in Table 1. As seen, 
the results indicated a satisfactory agreement with 
the literature reported solubilities [30—34]. 

The solubilities of uranyl nitrate at 60°C and 
two different pressures (150 and 250 atm) in pure 
and modified supercritical carbon dioxide were 
determined and the results are summarized in 
Table 2. Each reported solubility is the average 
value of three replicate samples. As expected, in 
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all cases, the solubility increases with increasing 
pressure [7—15]. 

Examination of the solubility data in Table 2 
reveals that, in the absence of any modifier, the 
solubility of uranyl nitrate in supercritical carbon 
dioxide at 60°C and 150 and 250 atm pressure is 
14.9 and 17.6 mg 1~', respectively. This observa- 
tion is in contrast to the case of toxic heavy metal 
ions whose solubility in pure supercritical CO, is 
reported to be negligible [35]. This is mainly due 
to relatively large ionic size of UO5* ion and its 
enhanced polarizability which in conjunction with 
NO, as a relatively hydrophobic anion will form 
a lipophilic pair which is easily soluble in an 
organic solvent such as methanol. This property 
has been the basis for many liquid—liquid extrac- 
tion methods for uranyl nitrate using organic 
solvent such as methylisobutyl ketone, diethyl 
ether and tri-n-butyl phosphate [36—40]. 

The data given in Table 2 reveal that addition 
of 50 ul methanol to the matrix (off-line method) 
will result in about 200—300 times enhancement in 
the solubility of uranyl nitrate in SF-CO,. It is 
well known that the main limitation of utilizing 
supercritical carbon dioxide in SFE is its inability 
to dissolve polar analytes even at very high densi- 
ties. However, solubilities of polar analytes in 
CO, can be greatly enhanced by increasing its 
solvating power via addition of various polar 
organic solvents known as modifiers. It is interest- 
ing to note that the direct delivering of the 
modifier to the sample matrix prior to the extrac- 
tion is shown to be more effective in promoting 
extraction efficiency than using in-line modifier 
addition [41,42]. Table 2 shows that addition of 
50 ul water to the sample matrix containing 
methanol will result in considerable decrease in 
solubility, compared to the case of addition of 
methanol alone. This is most probably because of 
the low solubility of water is SF-CO, [12], on one 
hand, and higher solubility of uranyl nitrate in 
water than that in methanol, on the other. 

From Table 2, it is seen that the addition of 50 
ul MeOH and 50 ul of | M nitric acid to the 
uranyl nitrate causes an exceptionally large in- 
crease in the solubility of uranium in SF-CO,, up 
to the values of 43.5 and 61.1 g1~!' at 150 and 
250 atm, respectively. Such enhanced extraction 


of uranyl ion in the presence of nitric acid with 
some organophosphorus reagents into SC-CO, 
has already been reported in the literature [15]. In 
nitric acid medium, not only the ion pair forma- 
tion of uranyl ion with nitrate anion tends to shift 
to the completion direction but also the solvated 
ion pair formation in supercritical fluid CO, will 
be facilitated. The lower solubility of uranyl ni- 
trate in the presence of acetic acid is most proba- 
bly due to reduced volatility of the uranyl acetate 
ion pair as it compared with uranyl nitrate. The 
observed reversed effect of hydrochloric acid on 
the solubility of uranyl ion in SF-CO, is due to 
the formation of charged complexed species be- 
tween the cation and chloride ion [43]. 


3.2. SFE studies 


To date, most of the published work on the 
successful SFE of metal ions has focused on the 
use of fluorinated ligands, which result in promis- 
ing recoveries of metals from liquid and solid 
materials [8—14,44]. In this work, the extraction 
of uranyl ion from cellulose matrix in the presence 
of a variety of chelating agents, belong to 
organophosphorous-, f-diketone-, macrocyclic- 
and amine-families, into methanol modified su- 
percritical CO, was investigated. The polar cellu- 
lose matrix used is not only simple and cheap, but 
also resembles other polar matrices in biological 
samples such as plants and animal tissues [11]. 
Moreover, because of both strong chemical and 
physical interactions between the solid matrix 
used and the solute, the extraction of ionic ana- 
lytes in SF-CO, could be a slow process. 

The extraction efficiency of pure and modified 
SF-CO, using different co-extractants for uranyl 
nitrate was first investigated and the results are 
given in Table 3. Without the co-extractant, some 
12% of the spiked uranyl nitrate could be ex- 
tracted from the filter paper by SF-CO,. This is 
significantly high, in comparison with the case of 
transition and heavy metal ions, the extraction of 
which in pure SF-CO, are reported to be neligible 
[10,11,35]. The result is well correlated with the 
significant solubility of uranyl nitrate in pure SF- 
CO, (Table 2). While, the addition of 50 pl of | 
M HNO; to the matrix will result in decreased 
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Table 3 
Percent extraction of uranyl ion from cellulose-based filter 
papers in the presence of 50 ul of different co-extractants* 





Modifier *% Extraction 





None 12.2 +2.5 
MeOH 25.6+ 0.4 
HNO, “5 

TBP 41.7412 
MIBK 34.7 + 2.0 
AA 26.2 + 2.2 





“Each filter paper is spiked with 100 pg U(VI) and 50 ul 
co-extractant. Conditions: T= 60°C and P = 250 atm. 


Table 4 

Percent extraction of uranyl ion from the cellulose-based filter 
papers in the presence of different chelating agents with 
methanol-modified SC-CO,* 








Ligand Percent extraction 





DBM IF 39 
HOX Seer eh 
DPASA 19.2+3.4 
DMPOPOP 14.3+ 2.6 
TOPO 49.0+1.8 
TOPO + DC18C6" 66.0 + 2.0 
HDEPH 63.2 + 3.2 
HDEPH + DC18C6° 70.8 + 2.5 
DC18C6 39.9 + 2.1 








“Each filter paper is spiked with 100 pg U(VI), 50 pg 
MeOH and 10 mg chelating agent. Conditions: T = 60°C and 
P = 250 atm. 

> Five milligram of each chelating agent. 


percentage extraction of uranyl nitrate to <5. 
This is most probably due to the strong interac- 
tion of acid with the cellulose matrix and the 
consequent enhancement in adsorption of uranyl 
ions on the solid support. However, addition of 
50 pul of other co-extractants will significantly 
increase the extraction efficiency of the system for 
uranyl ion in the order TBP > MIBK > AA > 
MeOH. The observed trend reflects the decreased 
complexing ability of the coextractants for UO3* 
ion from TBP to MeOH. Meanwhile, all of these 
compounds can also act as a polar modifier to 
increase the solvating power of CO, in the extrac- 
tion process. 

The results of extracting UO3+ from the cellu- 
lose-based filter papers in the presence of different 


chelating agents by methanol-modified SF-CO, at 
60°C and 250 atm are given in Table 4. It should 
be noted that, in all cases, the ligand was in large 
excess (10 mg) relative to the UO3* present in the 
sample. The organophosphorous reagents TOPO 
and, especially, HDEPH are known to not only 
form stable coordinatively solvated salts with 
uranyl nitrate through the P=O groups [15,45], 
but also possess sufficiently high solubilities in 
SC-CO, to be employed as a metal ion extractants 
in SFE [46]. Thus, these two ligands show the 
highest extraction efficiency for UO3* ion in the 
series. Other chelating agents used possess lower 
binding ability to UO3* ion and, therefore, result 
in decreased extraction of the cation into SF-CO,. 
On the other hand, the increased extraction effi- 
ciency of DBM and HOX over that of DPASA 
and DMPOPOP is partly related to the increased 
solubility of the former ligands in SF-CO, [47]. 
DC18C6 shows an extraction efficiency for UO3* 
ion which falls in between that of the organophos- 
phorous reagents and other chelating agents used. 

Synergistic extraction, which is frequently used 
to improve the extraction efficiencies of the lig- 
ands in conventional solvent—solvent extraction 
processes [32,48], is also shown to be applicable to 
the SFE systems [12,15,49]. The data given in 
Table 4 clearly illustrate the occurrence of such 
synergistic effect in SFE of uranyl ion. When 
TOPO and HDEPH are mixed with DC18C6, the 
extraction efficiency for uranyl ion significantly 
increased relative to that of the pure TOPO and 
HDEPH. As seen, the extraction efficiency for the 
TOPO + DC18C6 and HDEHP+DC18C6 sys- 
tems are 66.0 and 70.8%, respectively, exhibiting 
the highest percent extraction of uranyl ion into 
methanol-modified SF-CO,, among different lig- 
and systems used. 

Finally, the change in extraction efficiency of 
DC18C6 and HDEHP for uranyl ion in SF-CO, 
in the presence of 5% methanol as modifier by 
using the in-line method was studied and the 
results are shown in Table 5. As seen, in both 
cases, the percent extraction of uranyl ion shows a 
considerable increase, compared to that obtained 
by addition of 50 wl methanol via the off-line 
method. The recovery of uranyl ion under the new 
experimental conditions used by DC18C6 and 
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Table 5 ; 
Percent extraction of uranyl ion from the cellulose-based filter 


papers in the precence of DC18C6 and HDEPH with SF-CO, 
modified by 5% (w/w) methanol* 











Ligand Percent extraction 











DC18C6 +3 
HDEPH 97+4 











“Each filter paper is spiked with 100 pg U(VI) and 10 mg 
chelating agent. Conditions: T= 60°C and P = 250 atm. 


HDEPH jis 50 and 97%, respectively. A similar 
quantitative extraction of uranium from nitric 
acid solution [15] and from filter papers [12] into 
SF-CO, by the respective use of different 
organophosphorous reagents and 1:1 mixtures of 
TBP and fluorinated B-diketones has already been 
reported. In this case, it is very possible that 
interactions between the filter paper and the 
uranyl ion can be influenced by interactions with 
the ligands. Especially, it could be the reason that 
the in-line method is more efficient than the off- 
line method, where the amount of methanol to 
block the absorbing properties of the paper is 
considerably higher. 


4. Conclusions 


This study has demonstrated that the solubility 
of uranyl nitrate in SF-CO, in the presence of 
methanol and nitric acid is exceptionally high 
(some 61 g 1~'). The extraction ability of uranyl 
nitrate from cellulose-based filter papers into SF- 
CO, modified with TBP, MIBK, AA and MeOH 
was tested: the highest percentage extraction was 
obtained with TBP (some 42%). Percentage ex- 
traction of uranyl ion in the presence of eight 
different chelating agents HDEPH, TOPO, 
DC18C6, DBM, HOX, DPASA, DMPOPOP and 
TBP with methanol-modified SC-CO, was also 
investigated. The use of HDEPH resulted in the 
highest extraction efficiency (63%) in the series. 
Use of a 1:1 (w/w) mixture of HDEPH and 
TOPO with DC18C6 showed a synergistic effect 
in extraction of uranium into the supercritical 
fluid. Finally, the extraction of uranium into SF- 
CO, in the presence of 5% methanol by using the 


inline method was found to be quantitative with 
HDEPH as the most suitable chelating agent. 
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Abstract 








Polymerization of styrene in tetrahydrofuran (THF) and cyclohexane with CO, as antisolvent was carried out at 
335 K up to 10 MPa. The results showed that, with CO, as antisolvent, the molecular weights and molecular weight 


distributions of the product could be controlled by changing the pressure of the antisolvent. The polydispersity index 
of the polystyrene was smaller in comparison with the product polymerized in ordinary solutions. © 2001 Elsevier 


Science B.V. All rights reserved. 
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1. Introduction 


Supercritical (SC) CO, has strong ability to 
dissolve many non-volatile materials. Chemical 
reactions in SC CO, have been studied extensively 
[1-10], including polymerizations [1—6]. It is also 
well known that the solubility of many compounds 
in SC CO, is extremely low, i.e. it is a poor solvent 
for many substances. On the other hand, com- 
pressed CO, and some other gases are quite soluble 
in a number of organic solvents and expand them 
largely, which can reduce the solvent power of the 
solvents to such an extent that a solute dissolved 
is precipitated. This is the basis for gaseous anti-sol- 


* Corresponding author. Tel.: + 86-10-62562821; fax: + 86- 
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vent (GAS) process [11,12]. GAS processes have 
been successfully used in recrystallization of or- 
ganic solids [13,14], fractionation of natural prod- 
ucts [15] and preparation of ultrafine particles 
[16—22]. These kinds of processes have some unique 
advantages. For example, morphology, particle size 
and particle size distribution of the products can be 
tuned by pressure; dry particles can be produced in 
a single processing step; solvents can be recovered 
easily. 

The properties of some reactions vary with the 
properties of the solvents in which the reactions 
take place. The composition of a gas-liquid solution 
changes with the pressure of the gas, i.e. the 
properties of the solvent can be tuned continuously 
by pressure. In this work, styrene was polymerized 
in tetrahydrofuran (THF) or cyclohexane with CO, 
as anti-solvent. We aim at studying the 
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Fig. 1. Schematic diagram of the polymerization apparatus. l. 
Gas cylinder, 2,4,6. Valves, 3. High-pressure syringe pump, 5. 
Pressure gauge, 7. Temperature controller, 8. Optical cell, 9. 
Water bath, 10. Magnetic stirrer. 


effect of anti-solvent CO, on the properties of the 
products and partially replacing the organic sol- 
vents using environmentally acceptable CO,. 


2. Experimental 


2.1. Materials 


CO, with a purity of 99.9% was supplied by 
Beijing Analytical Instrument Factory. THF, sty- 
rene, cyclohexane and 2,2’-azobis(isobutyronitrile) 
(AIBN), supplied by Beijing Chemical Factory, 
were A.R. grade. Styrene was washed twice with 
10% NaOH solution and finally with water until it 
was free from alkali, dried over fused CaCl, for 24 
h and distilled under reduced pressure. AIBN was 
recrystallized twice with methanol and stored in the 


dark in a refrigerator. 


Table | 





| without CO, 





@ with CO, as antisolvent 








10 20 30 40 


50 


Initial Monomer Concentration / % 


Fig. 2. The yields of PS polymerized with and without CO, as 


antisolvent. 
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Fig. 3. The Mw of PS polymerized with and without CO, as 


antisolvent. 


2.2. Polymerization 


The apparatus is shown schematically in Fig. 1. 
It was composed mainly of a gas cylinder, a 
high-pressure syringe pump, an optical cell, a 
pressure gauge and a constant temperature water 
bath. There were graduations on the optical cell, 


Mw, Mn, PDI and yields of PS polymerized in THF expanded with different pressure of CO,* 





C., (%) 


P (MPa) 


VER 


Mw x 10~4 (g.mol~') 


Mn x 10~4 (g.mol~!) PDI 


Yields (%) 








50 
33 
25 
17 
10 


0.0 
5.8 
7.8* 
9.1% 
10.0* 


1.0 
1.5 
2.0 
3.0 
5.0 


8.08 
4.21 
2.26 
1.89 
1.27 


3.63 
2.30 
1.44 
1.42 
1.00 


2.2 
1.8 
1.6 
1.3 
1.3 


59 
41 
35 
24 
1] 








“Reaction condition: 335 K, 16 h, the concentration of AIBN (after expansion with CO.) is 0.06 mol.I~'. C.,, is the initial 


monomer concentration in volume percentage after CO, expansion; VER, volume expansion ratio; Mw, weight-average molecular 
weight; Mn, number-average molecular weight; PDI, polydispersity index of the molecular weight distribution; Yields were 


determined gravimetrically. aa 
* PS precipitated out when the reaction time is long enough. 
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Fig. 4. The Mn of PS polymerized with and without CO, as 
antisolvent. 
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Fig. 5. The polydispersity of PS polymerized with and without 
CO, as antisolvent. 


so that the volume of the liquid in the cell could 
be known easily. 

In a typical experiment, suitable amounts of 
initiator (AIBN), liquid solvent, and monomer 
were charged into the optical cell at room 
temperature. The cell was then transferred into 
the water bath at 335 K. After thermal 
equilibrium was reached, CO, was charged into 
the cell and magnetic stirrer was started. The 
volume of the liquid was expanded continuously 


Table 2 
Mw, Mn, PDI and yields of PS polymerized in THF* 


because CO, was dissolved. This was continued 
until the desired liquid volume was _ reached, 
which could be easily known from the 
graduations on the optical cell. The reaction 
lasted for 16 h. At the end of the reaction, the cell 
was cooled, CO, in the cell was slowly vented and 
the polymer product was removed from the cell. 
After the venting of CO,, the mixture was 
homogeneous. At some conditions, the viscosity 
of the mixture was very high, so it was rinsed with 
THF. The polymer was precipitated in methanol, 
separated and then dried in vacuum. 


2.3. Characterization 


The weight-average molecular weight (Mw), 
number-average molecular weight (Mn) and the 
polydispersity index of the molecular weight distri- 
bution (PDI = Mw/Mn) were determined by gel 
permeation chromatography (GPC) using UV-Vis 
detector (THF as solvent; flow rate, 1 ml/min; 
injection volume 10 ul; detecting wavelength, 254 
nm). The yield was determined gravimetrically. 


3. Results and discussion 
3.1. Fixed monomer /solvent ratio 


Styrene was polymerized in THF with and with- 
out CO, anti-solvent. The concentration of initia- 
tor AIBN was controlled to be 0.06 mol/l after the 
expansion with CO,. The volume ratio of THF/sty- 
rene was fixed at 1:1 and the pressure of CO, was 
varied from 0 to 10 MPa. The volume of the liquid 
was expanded by the dissolved CO,. The volume 
expansion rate (VER), 1.e. the volume ratio 








Cin (%) Mw/10* (g.mol~') 


Mn/10* (g.mol~ ') PDI 


Yields (%) 








50 8.08 3.63 
33 5.39 2.62 
16 1.68 1.23 


an0 59 
2.0 47 
1.4 34 








4 Reaction condition: 335 K, 16 h, the concentration of AIBN (after expansion with CO,) is 0.06 mol.l~'. C,, was the initial 
monomer concentration in volume percentage after CO, expansion; Mw, weight-average molecular weight; Mn, number-average 
molecular weight; PDI, polydispersity index of the molecular weight distribution. Yields were determined gravimetrically. 
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Table 3 
Results of styrene polymerization in CO,-THF (the initial volume concentration of the monomer in the expanded solution was 


33%)* 





P (MPa) VER Mw/10* (g.mol~') Mn/10* (g.mol~') PDI Yields (%) 





0.0 1.0 be 1.74 59 
3.1 1.2 3.28 1.67 5] 
5.8 1.5 ear 1.44 48 








* Reaction condition: 335 K, 16 h, the concentration of AIBN (after expansion with CO,) is 0.06 mol.l~'. VER, volume expansion 
ratio; Mw, weight-average molecular weight; Mn, number-average molecular weight; PDI, polydispersity index of the molecular 
weight distribution; Yields were determined gravimetrically. 


Table 4 
Results of styrene polymerization in CO,-cyclohexane (the initial volume concentrations of the monomer in the expanded was 


solution 50%)* 








P (MPa) VER Mw/10* (g.mol~') Mn/10* (g.mol~') Mw/Mn Yields (%) 








0.0 1.0 21.35 9.19 2.32 60 
6.0 1.25 10.13 4.59 2.21 47 
8.4* 1.67 7.74 3.84 2.02 4] 
9.4* 2.0 6.57 3.34 1.96 34 








* Reaction condition: 335 K, 16 h, the concentration of AIBN (after expansion with CO,) is 0.06 mol.1~'. VER, volume expansion 
ratio; Mw, weight-average molecular weight; Mn, number-average molecular weight; PDI, polydispersity index of the molecular 
weight distribution; Yields were determined gravimetrically. 

* PS precipitated out when the reaction time is long enough. 
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Fig. 6. Effect of VER on the yield, molecular weight and its distribution of PS polymerized in THF-CO, and cyclohexane-CO, 
(CH-CO,). 
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of the liquid after and before dissolving CO,, was 
a function of pressure. The initial volume con- 
centration of the monomer (concentration before 
reaction) after the expansion with CO, (C,,) 
changed from 50 to 10%. 

The results of the yield, Mw, Mn, and PDI of 
polymerization at different VER are shown in 
Table 1 and Figs. 2—5. It can be seen that the Mw, 
Mn, PDI and yield decrease with increasing 
pressure of CO, or VER, i.e. the molecular weight 
and molecular weight distribution can be con- 
trolled by pressure. There are two possible reasons 
for this. Firstly, the property of the solvent changes 
with VER. Secondly, the initial monomer con- 
centration (concentration before reaction) in the 
expanded solution decreases as VER is increased. 
In order to study which factor is dominant, 
polymerization was conducted in pure THF with 
different initial monomer concentrations. The 
solution can be considered as expanded with THF 
as initial monomer concentration is reduced. The 
dependence of average molecular weight, PDI and 
yields on VER is shown in Table 2 and Figs. 2—5. 

It can be seen from Figs. 2—5 that the Mw, Mn, 
PDI and yields of the polystyrene (PS) polymerized 
in pure THF also decrease with decreasing initial 
monomer concentration. However, the Mw, Mn, 
PDI and yields of the product polymerized in 
CO,-THF are lower than those polymerized in 
THF when the monomer concentration is the same. 
This indicates that the decreasing of the yields, the 
molecular weight and molecular weight distri- 
bution of PS with the addition of CO, result from 
not only the effect of the dilution of monomer, but 
also from changing the properties of the 
polymerization medium. 


3.2. Fixed initial monomer concentration 


In the above section we discussed the effect of 
volume expansion on the polymerization at fixed 
monomer/solvent ratio. We also conducted the 
polymerization in CO,-THF mixtures with differ- 
ent styrene/THF ratios while the initial monomer 
volume concentration after expansion was con- 
trolled to be 33%. The concentration of AIBN 
was also 0.06 mol/l after the expansion with CO,. 
The results are shown in Table 3. The polymeriza- 


tion was also carried out in CO,-cyclohaxene with 
different styrene/cyclohexane ratios and the initial 
monomer volume concentration was 50%. The 
results are given in Table 4. 

The data in Tables 3 and 4 indicate that Mw, 
Mn, PDI and yields also decrease with increasing 
VER or pressure of CO, at fixed monomer con- 
centration. These results support the conclusion 
that the decreasing of the Mw, Mn, PDI and 
yields of PS with addition of CO, is not only the 
result of the dilution of monomer. The properties 
of the product can be tuned by the pressure of 
anti-solvent CQO,. 


3.3. Effect of the anti-solvent CO, 


All the polymerizations started in homogeneous 
solutions and the solutions remained to be homo- 
geneous with reaction time at the lower pressures. 
At higher pressures, however, the polymerization 
changed from homogeneous to heterogeneous as 
the polymerization proceeded. There are mainly 
two reasons for this. The first is that CO, is a 
poor solvent and the concentration of CO, in the 
solution increases with pressure. Secondly, the 
monomer is a good solvent for the polymer and 
its concentration is reduced with reaction time. In 
Tables 1 and 4 we marked the pressure at which 
precipitation of PS occurs when the reaction pro- 
ceeded long enough. 

In organic solvents the gel effect may appear in 
solution polymerization as the solvent changes 
from good solvent to poor solvent. It results in 
increasing of the molecular weight of the prod- 
ucts. The gel effect is the result of increased 
viscosity, which reduces the rate of the terminal of 
the radicals. In our experiments, the solvent 
power decreases as the VER increases. However, 
gel effect was not found, although the solvent 
powers of the solvents were reduced by dissolved 
CO,. The possible reason is that the viscosity of 
the solution with expanded with CO, is low. 


3.4. The effect of the solvent power 


It can be seen from Tables 3 and 4 that the 
molecular weight of the PS polymerized in THF 1s 
far less than that of polymerized in cyclohexane. 
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However, the yields are almost identical for the two 
solvents. It may partially attribute to the difference 
in initial monomer concentrations. On the other 
hand, the different solvent power of THF and 
cyclohexane also result in the different properties 
of the PS. THF is a good solvent for the PS while 
cyclohexane is a relatively poor solvent. 

The variation of yields, molecular weight and 
molecular distribution with VER in the two solu- 
tions is shown in Fig. 6. This figure indicates that 
in both solvents the molecular weight and its 
distribution decrease with increasing VER. 


4. Conclusion 


In this paper, we have conducted the polymeriza- 
tion of styrene in THF and cyclohexane expanded 
with CO,. The results showed that the yields, 
molecular weight, molecular weight distribution 
and polydispersity index of the polystyrene could 
be tuned by controlling the pressure of CQ . 
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